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[1] Mathematikunterricht mit Computeralgebra-Systemen, H. Heugl, W. Klinger und J. Lechner 
307 Seiten, DM 59,90; oS 443.00; sFr 48,00; Addison- Wesley, 1996, ISBN 3 8273 1082 2 

Ein didaktisches Lehrerbuch mit Erfahrungen aus dem osterreichischen DERIVE-Projekt. Eine 
Kurzbesprechung folgt im nachsten DNL. 

[2] TI 92 - du lycee a la prepa, Henri Lemberg 

DUNOD / Texas Instruments France, Paris Mai 1996, ISBN 2 10 003039-6 

This book recapitulates on 314 pages the mathematics curriculum of French gymnasiums and 
shows up all the concepts and algorithms which a French student needs to know for the entrance 
examination for French universities. 

[3] TI 92 - les programmes!, Jean-Michel Ferrard 

DUNOD / Texas Instruments France, Paris Mai 1996, ISBN 2 10 003104-X 

This is a real treasury of 480 pages full with TI 92 functions and programs. All the programs can 
be used for their own but they also can be assembled within a library divided in different folders 
containing their own menus which allow an easy approach to the programs. The sections co- 
vered are: Arithmetique et Trigonometrie, Poiynomes, Matrices, Geometrie, Developpements Li- 
mites, Analyse, Geometrie Differentielle, Fonctions Speciales, Probabilites. Even if you only have 
a very poor knowledge of French - like me - it is easy to follow. Fortunately the language of maths 
is international. If you want to learn and to train programming with the TI 92 then I'd strongly rec- 
ommend this book. 

[4] TI 92 - le "top" des jeux!, Vincent Bastid et Emmanuel Neuville, DUNOD / Texas Instruments France, 
Paris Mai 1996, ISBN 2 10 003040-X 

Transformez votre calculatrice TI 92 en console de jeux! 

A collection of games containing Tl-tris, Demineur (Minesweeper), Bataille-Navale (for two TIs!), 
Tl-lnvaders, Tic-Tac-Toe-3D and four more games. 

I have to thank Mathias Makowsky from Marbach, Germany, who faxed the titles of the three 
French books. Fie saw them in a bookstore in Brittany during his holidays. Many thanks, the 
books you recommended are very useful. Josef 

[5] An Introduction to the Mathematics of Biology (With Computer Algebra Models) 

Yeagers, Shonkwiler and Herod, Birkhauser Boston 1996, ISBN 0 8176 3809-1 
You can find a short review on page 1 7. 

[6] AGNESI to ZENO, Over 100 Vignettes from the History of Math, Sanderson Smith 
Key Curriculum Press, Berkeley 1996, ISBN 1 55953 107 X 

This is not a Computer Algebra Book, but it is a wonderful book to motivate teachers and stu- 
dents as well for investigations and projects and presents a lot of facts concerning history of 
mathematics and the men and women who wrote this history. The book is available from Jan 
Vermeylen's Rhombus Shop. (See the address on page 34). 


Call for partners 

At the Information Day about the European Union’s Educational Programs I was 
asked eventually to start a transnational project in the frame of the COMENIUS Pro- 
gram. So I’d like to ask for partners who could imagine cooperating. We have to be 
at least three participants from different European Union countries. I have two ideas 
in mind (but maybe there are better ones): 

• Set up a transnational structure for math teacher's training in modern technolo- 
gies. 

• Exchange, evaluate, improve and customize teaching materials for teaching maths 
using modern technologies. 

If you are interested, then please contact me as soon as possible. Please notice my 
new email address: Josef Boehm@bboard.at (is not valid since a couple of years!) 


LETTER 


THE EDITOR 


Liebe DUG Mitglieder, 

Einige recht arbeitsreiche Tage und Ndchte liegen hinter 
mir. Aber es ist wieder gelungen: der letzte DNL des Jahres 

1996 ist fertig. Wahrend er nun zum Drucken geht, wird 
noch rasch die Diskette des Jahres 96 randvoll gepackt und 
iiberpriift, ob auch alle Dateien und „ Weihnachtsgeschen- 
ke “ drauf sind. Dann werden meine Frau Noor und dieses 
Mai auch meine Tochter Astrid einige Hundert Kopien 
ziehen, die Newsletter mit Diskette, 3D-Brille und Renewal 
Form in ein grofien Kuvert stecken, Uberall die notwendi- 
gen Stempel anbringen und sie dann auf die weite Reise 
schicken. 

Bitte beachten Sie meine neue email- 
Adresse am Ende dieser Seite. Endlich 
habe ich meinen eigenen Internetzugang 
an der Schule. Dass dieses Medium bereits 
geniitzt wird, zeigt das reiche User Forum. k 

Es bietet sich auch ein neues Service der 
DUG an: Falls Sie den einen oder anderen 1 
Artikel aus einem DNL - auch von friihe- 
ren - als Textfile brauchen konnten, kann 
ich Ihnen diesen gerne iiber email schi- 
cken. Heute habe ich wieder in meinem 
elektronischen Postkasten gefischt und 
hatte Anglergliick : es gibt eine Antwort 
von Al Rich auf das Matrizenproblem im 
User Forum und Terence Etchells machte eine aufregende 
Ankiindigung fur einen moglichen Beitrag im nachsten 
Jahr. 

Ich mochte Sie auch nochmals auf meine Partnersuche fiir 
ein EU-Projekt im Rahmen des COMENIUS Programms 
auf der Infoseite aufmerksam machen. Ich konnte heuer 
zwei Klassen mit TI 92 ausstatten, daher ware ich an einem 
Austausch an Unterrichtsmaterialien fiir den TI 92 - aber 
nach wie vor auch fiir DERIVE - sehr interessiert. Ich 
arbeite heuer besonders im Precalculus Bereich 

Meinen letzten Brief dieses Jahres mochte ich aber nicht 
beenden, ohne Ihnen alien fiir die engagierte Mitarbeit zu 
danke. Ohne diese gabe es Icing st keinen DNL me hr. Noch 
wertvoller aber finden wir die herzliche Freundschaft, die 
die DERIVIANER weltumspannend untereinander verbin- 
det. Jedes Zusammentreffen, jede Konferenz ist ein deutli- 
ches Zeichen dafiir. Viele personliche Begegnungen und 
Freunds chaf ten haben sicher nicht nur unser (Noor und 
Josef) Leben, sondern das von vielen unter uns bereichert. 

In diesem Sinne wiinschen ich alien ein frohes Weihnachts- 
fest und ein gesundes, gliickliches und erfolgreiches Jahr 

1997 - Jahr 1 nach DERIVE for Windows. 

Ich hoffe Euch alle im nachsten Jahr ^ 
wieder begriifien zu diirfen. ( 



Dear DUG Members , 

Some very busy days and nights are lying behind me. 
But it has lucked again: the last DNL of 96 is ready. 
While it is going to be printed, the Diskette of the Year 
has to be filled brimful and checked if all the files and 
the Christmas gifts are on it. Then my wife Noor and 
my daughter Astrid will produce several hundreds of 
copies, put them together with the DNL, the 3D- 
spectacles and the renewal form in big envelops, affix 
the stamps and address labels and send them on their 
big journeys. 

Please note my new email address at 
the end of the page. I am happy to 
have my own email access now at my 
^ school. And I like to use this media, as 

Jj\ you can see in the rich User Forum. So 

jR Td like to offer another DUG service: if 

\ you would like to have the text file of 

* any DNL article - even from earlier 

|) DNLs - then I could email it to you. 

Today I was fishing again in my elec- 
tronic mailbox, and I was lucky. Among 
^ others I found Al Rich's answer to the 
matrix problem in the User Forum and 
my friend Terence Etchells made an 

exciting announcement for a possible 

contribution for 1997. 

I want to focus your attention once more on my call for 
partners for an EU-project within the COMENIUS 
program on the Information page. I was able to equip 
two of my classes with TI 92s and so I am very inter- 
ested to exchange teaching materials for the TI - but 
also for DERIVE. In these courses I deal mainly with 
precalculus stuff. 

I don't want to finish my last letter of 96 without thank- 
ing you all for your enthusiastic cooperation. Without 
this there would be no DNL any longer. But we find 
much more valuable the warm hearted friendship 
which connects the DERIVIANS from all over the 
world. Each meeting, each conference is a witness of 
this fact. Many personal friendships and meetings 
have not only enriched our lives (Josef and Noor) but - 
I am sure - those of many of us. 

In that sense I wish you all a Merry Christmas and a 
healthy, happy and successful 1997 - year 1 after 
Derive for Windows. Hope to meet you all next year 
again! 

Sincerely yours ^ r 


(Please settle your membership fee for 1997) 


(Bitte begleichen Sie Ihren Mitgliedsbeitrag 1997) 


email: nojo.boehm@pgv.at 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 40 pages minimum. The goals of 
the DNL are to enable the exchange of ex- 
perience made with DERIVE as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

We include now a section dealing with the 
use of the TI-92. 


Editor: Mag. Josef Bohm. 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone. 43-(0)660 3136365 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
D-N-L. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Preview: Contributions for the next issues 

LOGO in DERIVE, Lechner & Roanes Lozano, AUT & ESP 
3D-Geometry, Reichel, AUT 
Algebra at A-Level, Goldstein, UK 

Graphic Integration, Linear Programming, Various Projections, Bohm, AUT 

Tilgung fremderregter Schwingungen, Klingen, GER 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Linear Mappings and Computer Graphics, Kiimmel, GER 

Solving Word problems (Textaufgaben) with DERIVE, Bohm, AUT 

Line Searching with DERIVE, Collie, UK 

About the "Cesaro Glove-Osculant", Halprin, AUS 

Tangrams with DERIVE, Poblacion, ESP 

Hidden Lines, Weller, GER 

Fractals and other Graphics, Koth, AUS 

Experimenting with GRAM-SCHMIDT, Schonefeld, USA 

The TI-92 Section, Waits a.o. 
and 

Setif, FRA; Vermeylen, Belgium; Leinbach, USA; Halprin, AUS; Biryukow, RUS; 
Weth, GER; Wiesenbauer, AUT; Keunecke, GER; Aue, GER; 

Stahl, USA; Mitic, UK; Sirota, RUS;and 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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DNL: Some days later there was another message from Terence which might be interesting 
for pure DERIVIANS and Tl-Users as well: 

How are you? It was nice to meet up with you again in Bonn. (I enjoyed our talks , too) 

Whilst discussing the TI 92 with David Stoutemyer at Bonn I was bemoaning the fact that the Ti did 
not have my favourite function ITERATES. He said ”no problem, we’ll write one”: Dave then gave 
me a brief tutorial on writing functions in the TI 92 function programming language. Time was short 
so he scribbled a few ideas on paper. On the trip back from Bonn and over Summer I set to write a 
series of DERIVE functions for the TI 92, such as: ITERATES*}, ITERATE*}, ELEMENT*}, 
DEL_ELEM(), (This is DELETEELEMENTQ, but the TI 92 restricts function names to 8 charac- 
ters), SWP_ELEM(), REV_ELEM(), RHS(), LHS(). Also, DERIVE will easily plot a2x/i matrix in a 
2D-plot window. I could not find an easy way to do this with the TI 92 so I wrote a program 
PLOTMAT that plots a2x« matrix. 

I will, very shortly, be putting these functions on my Web page; see the signature file below for the 
URL if you should wish to browse it (there are lots of my files to down load as well). 

e mail: t.a.etchells@livjm.ac.uk 

web page: http://www.cms.livjm.ac.uk/www/homepage/cmstetch/index.htm 

DNL: It would take a lot of space to print all the functions. You will find them on the Diskette 
of the Year in subdirectory TI, accompanied by W.Propper's DIRA package for investigating 
functions and his SOLSYST-function to solve simultaneous equations in a DERIVE like 
way. That are the true DERIVIANS, who implement their favourite DERIVE functions on the 
TI-92. (See more in the TI-92 Corner and on page 43- a demo of PLOT_MAT.) 


Johann Wiesenbauer, Vienna, Austria 

In the following the polynomial system of equations: 

(cf. The International DERIVE Journal, Vol.3, No. 2, p. 96) 
will be solved by means of my routines (cf. DNL#23, Titbits(8)). 

(Johann refers to an article dealing with the implementation of Groebner bases in DERIVE. 
Josef) Preload RED(u,v) and S0LVE2(u,v,x,y) from TITBITS8.MTH 

First we prove that there doesn’t exist a solution of the system above where x, y, z are all different. 

The following equation 

2 

#1: RED(x*z - 1 + y, y-z -1 + x) = y (1 - z ) + z - 1 

shows that either z = 1 or y = l/(z + 1). For reasons of symmetry z = 1 or x = l/(z+l) holds as well. 
Therefore our assertion is certainly true for z ^ 1 . But for z = 1 the system above simplifies to x y = 0, 
x + y = 1, which leads to x = z or y z due to 



#2: S0LVE2 (x • y , x + y = 1) 



U ser=Simp (User) 


Assuming w.l.o.g. thaty = z we are left with 


2 

#3: S0LVE2(x*y - 1 I y, y -1 + x) 


0 

V5 1 

2 2 

V5 1 


1 

V5 1 

2 2 

V5 1 


2 


2 


2 


2 
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The other two cases x=y and x = z, respectively, are obtained by simple permutations of the solutions 
above. They yield two additional solutions, namely [1, 0, 1] and [1, 1, 0], TITBITS 9 deals also with 
this problem. 

SOLUTIONS of DERIVE 6 and solve of the TI-92 have no problems solving this non linear 
system, Josef 2010 


SGLUTIQNS(x-y = 1 - z a x ■ z = 1 - y a y ■ z = 1 - x , [x , y , z]) 


Oil 
10 1 


110 


A 1 Is 1 IS 1 

2 2 2 2 2 2 


Is 1 Is 1 IS 1 

2 2 2 2 2 2 . 


Fi7»mjY FF? Y f:=t t fht y ft Y Fe^ Y 'i 
t f — |fl 1 gebra |Ca 1 c |oiher |Pr gn 1 0 |c 1 ean Up | | 


1 solve(x ■ y = 1 
-Gf5 + l) 


■ z and x ■ z = 1 


and y = 


y and y ■ z ► 

(J5 - 3) (J5 + l) 3 ► 

16 anri 


...d b=1 01 * x=Q and <j=l and z=l 

MAIN RAD EXACT FUNC 1/30 


Heinz Rainer Geyer, St.Katharinen, Germany 

Vielen Dank fur Deine Materialien. Bei meinen Bemiihungen um „Teilermengen“ habe ich eine Funk- 
tion SET und die entsprechenden Mengenoperationen in DERIVE vermisst. Umstandlich war auch, 
dass die IF-Funktion keine echte Leerausfuhrung zulasst. Das ”?” ist schwer zu bearbeiten 

Moglicherweise lasst sich die direkte Farbwahl ahnlich wie in BASIC direkt implementieren. 

In meiner 6. Klasse bin ich gerade bei den Teilem und Vielfachen. Ich habe etwas herumexperimen- 
tiert. Hier findest Du die Ergebnisse meiner Zusammenfassung: 


#1: 

# 2 : 

#3: 

#4: 

#5: 

#6: 

#7: 

#8: 


Set of divisors, 1st attempt: — > Jn 


1 

r 

" n | 


1 1 

> 


DIVISOR l(n) :z VECTOR 

1 

l iF 

r i 

= FEOORj 

f n 1 

It] 

f n 1 

|, i, 0 
1 1 

, i, 1, Jn 

J 


DIVISOR2(n) :r VECTOR 

H 

1 = FLOOR | 

It] 

|. i, oj 

, i, FLOOR (,/n) + 1, n 

j 


PARTl(n) := SELECT(k t 0, k, DIVISORl(n)) 

PART2(n) := SELECT(k t 0, k, DIVIS0R2(n)) 

DIVSET(n) :i APPEND(PARTl(n) , PART2(n)) 

DIVSET(3450) 

[1, 2, 3, 5, 6, 10, 15, 23, 25, 30, 46, 50, 69, 75, 115, 138, 150, 230, 345, 575, 690, 1150, 
1725, 3450] 


#9: 

# 10 : 

# 11 : 


0.11 sec - in 1996 10.8 seconds only by approximating #7 
Attempt #2: Why not counting until n? 

DIVISORO(n) := VECTOR 


Ll_ 
1— 1 

11 1 
= FLOOR 


, i, 0 

. i i 1. n 

s. 

l i 1 

l i J 

1 J 

1 J 


# 12 : 


# 13 : 


PARTO(n) := SELECT(k t 0, k, DIVISORO(n)) 
PART0(3450) 
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#14: 


[1, 2, 3, 5, 6, 10, 15, 23, 25, 30, 46, 50, 69, 75, 115, 138, 150, 230, 345, 
1725, 3450] 


575, 690, 


1150, 


#15: Again 0.11 sec - in 1996 it needed 10.3 sec 

#16: Attempt #3: it should work faster if working with half of the divisors! 

f " "I 

#17: PART3(n) := VECTOR , k, PARTl(n) 

l k J 

#18: PART3(3450) = [3450, 1725, 1150, 690, 575, 345, 230, 150, 138, 115, 75, 69] 

#19: D_S ET ( n ) :z APPEND(PARTl(n) , R EVERS E( PART3 ( n ) ) ) 

#20: D_S ET(3450) = [1, 2, 3, 5, 6, 10, 15, 23, 25, 30, 46, 50, 69, 75, 115, 138, 150, 230, 345, 
575, 690, 1150, 1725, 3450] 


#21: seems to work, but ...: 

#22: D_SET(625) = [1, 5, 25, 25, 125, 625] 

#23: ... hence: 

DSET(n) := 

If -Jn t FlOOR(Vn) 

#24: D_SET(n) 

DEEETE_EEEMENT(D_SET(n) , DIM(D_SET(n))/2) 

#25: DSET(3450) = [1, 2, 3, 5, 6, 10, 15, 23, 25, 30, 46, 50, 69, 75, 115, 138, 150, 230, 345, 

575, 690, 1150, 1725, 3450] 


#26: DSET(625) = [1, 5, 25, 125, 625] 


Ich bin sicher, es gibt bessere Losungen. Jedenfalls ist mir der Unterschied zwischen Simplify und 
approX deutlich geworden. 

(Heinz Rainer complained that there are no set operations in DERIVE, he also would like to 
have a true "not - execution", the "?" is not very comfortable to work with. He also has the 
idea to choose the plot colour directly by a command (similar to BASIC). Working with divi- 
sors and multiples in his 6 th form he experimented with sets of divisors. See Rainer’s results. 
He is sure that there are better solutions, but he has learned the difference between Simplify 
and approX. In DERIVE 4.x you will find Set operations. And the DIVISORS are also imple- 
mented.) 


DIVISORS(3450) 


= [1, 2, 3, 5, 6, 10, 15, 23, 25, 30, 46, 50, 69, 75, 115, 138, 150, 230, 345, 
575, 690, 1150, 1725, 3450] 


George Freeman and A1 Kopf, Kuweit City 

Dear Josef, we have come across this bug in DERIVE 3.02: 


# 1 : 


F(x) 1 G(x) 

2 

4(1 + X ) 


x 


2 

4(1 - x ) . 


#2: G(F (x)) = x 


2 

#3: F(G (x)) = - x-SIGN(x - 1) 


Can you tell us why? 
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DNL: Try the following: 


#4: x :e Real (-1, 1) 

#5: G(F(x)) = x 

#6: F(G(x)) = x 

I believe I know George and Al. Kopf seems not to be an Arabic name. Yussuf. 


Rudiger Baumann, Celle, Germany 

... Mit ihrer Punktfolge haben Carl Leinbach & Marvin Brubaker (DNL#22 rev, page 29) ein hiibsches 
Beispiel geliefert. Verallgemeinert man die Konstruktion ein wenig, gelangt man zu Edward Sawada's 
"misguided missile" (DNL#22 rev, page 8): 

"MISMIS2 . MTH" 

[pO : = [ 0 , 0 ] ,pl:=[0.5,0. 5*SQRT ( 3 ) ] ,p2:=[l,0] ] 

P (r, s, k) : =IF ( k=0 , pO , IF ( k=l , pi , IF ( k=2 , p2 , r*P (r, s, k-3) +s*P (r, s, k-2) ) ) ) 
FOLGE1 (r, s, n) :=VECTOR (P (r, s, k) , k, 0, n) 

FOLGE1 (0.9,0.1,15) 


Wegen der entsetzlich langsamen Rekursion empfiehlt sich eine iterative Fassung: 

FOLGE2 (r , s, n) : =ITERATES ([b,c,r*a+s*b], [a, b, c] , [p0,pl,p2],n) 

Der Aufruf von 
FOLGE2 (0.9,0.1,15) 

liefert das ominose Sawad'sche "missile". Interessant sind auch die Falle r + s > 1 und r + s < 1, bei- 
spielsweise folge2(0. 15,0.9,50) und folge2(0.05,0.9,80). Hier konnen die Schuler selbstandig experi- 
mentieren und Vermutungen hinsichtlich Konvergenz bzw. Divergenz aufstellen. 

Rudiger Baumann shows a link between Carl Leinbach’s & Marvin .Brubaker’s recursive se- 
quence of points in DNL#22 rev, p. 29 and E. Sawada’s "misguided missile" on p. 8 (same 
issue). Generalizing the construction presented in the first contribution we obtain the "mis- 
sile". As the recursive construction is very slow, Rudiger recommends the iterative proce- 
dure. Furthermore he points out that there are interesting cases with r+s>1 and r+s<1. 
Students could set up conjectures concerning convergence and divergence. 


Alfonso J. Poblacion, Valladolid, Spain 

Hello dear DERIVERS, 

Try to plot the function x A (l/2) (I did using DERIVE 2.56 and 3.0). You can see that, obviously, this 
function takes [0,+inf) — > [0,+inf). Now plot ABS(x A (l/2)), and surprisingly, you will have a function 
defined in the whole real line. How can it be? 

Then if you Simplify ABS(x A (l/2)), DERIVE gives you ABS(x A (l/2)). So, if you Simplify 
ABS((-4) A (l/2)), the result is 2 and this is not right. And the same happens for every even root. 

I think this happens because DERIVE has implemented this simplification for odd roots for which it is 
valid. But I think this must be modified for even ones, don’t you? 

If I did something wrong or anyone has a further explanation, please let me know. Thanks. 
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DNL : There were some explanations in DERIVE News Group. In general all had the same content, so 
I'll combine the different answers (michel-gosse@magic.fr, Al Rich, A. van der Meer, Josef, a.o): 

The problem with DERIVE is that it calculates with complex numbers. The absolute value of 2i then 
becomes 2. 

Like virtually all CAS, DERIVE works in the complex domain, not just in the real domain. 

You can define a function that will plot correctly, but that is a little tricky: 

f(x) := if (x < 0, 0, abs (x A (1/2) ) ) . 


Alfonso J. Poblacion, Valladolid, Spain 

Thank you very much about your message 
explaining what DERIVE does with the ABS- 
function. All the answers pointed out in the 
same direction, including one that I discovered 
later in the DNL# 14, page 6, messages 2708 
and 2710. I knew (but I forgot) that DERIVE 
works with complex numbers (that was a terri- 
ble oblivion), but what I did not know was that 
it plots in the way it does with these functions. 
A lot of our students usually work with 
DERIVE alone, and they can be confused 
about these behaviours, in case they detect 
them. Imagine a complicated function that 
involves ABS and SQRT and its plot: they 
believe in what they are seeing, so we must 
advise them (also, most of them do not know 
so much about complex numbers. In their first 
course they only deal with real numbers). 


I have got a lot of tangrams. I will try to send 
them by e-mail 



DNL: Thanks for the tangrams. I'm sure they will show together with your contribution a new 
and still unknown facette of DERIVE and they combine in an ideal way mathematic tuition 
with entertainment. 


Sergey Biryukow, Moscow, Russia 

Dear Josef, I have written functions for Implicit 3D 
Plots in DERIVE and I am writing DOC & DMO 
files now. Functions for Implicit 3D Plots in 
ACROSPIN are also written and return a vector of 
lines (2x2 numeric matrices). Is this format 
compatible with the tool you are going to present in 
the next DNL? I shall try to polish my IMPLICIT 
PLOTS utility as fast as possible in order to include 
it in the 1996 DNL diskette. Sincerely Sergey. 

x 2 y 1 + y 2 z 2 + x 1 z 2 = 1 
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Harald Lang, Stockholm, Sweden (lang@math.kth.se) 

My version of DERIVE Classic 3.13 has a bug. If I try to solve 


(z - #i)/(z + #0 - (w - 1 )/(w + #0 = 0 


for w, it says "memory full" after a few seconds. However, solving for z is fine. Both z and w are de- 
clared Complex. Cheers — Harald Lang 


Answer from Alain Pomirol , Langon, France (Pomirol@aol.com) 

Same problem with 3.1 1 version of DERIVE XM. The given equation is simplified to 

z/(z + #/) + (1 -w)l(w+ #/) - #//(z + #/) 

but is not factorized. The solution is: 

1 . Factorize 

2. Solve (equation, w) 

It's a problem of the employed method. Note: The Tl 92 csolves this equation without any 
problems. Yours truly - Alain Pomirol 

DNL : Same problem in DfW. Josef 


Z - L W - L 


#1: 

Z + L 1 

#2: 

z :€ Complex 

#3: 

w :€ Complex 

#4: 

SOLVeI 

1 z + 


( Z - i W-L ^ 

=0 ( W = 

V Z + L W + L ) 


Solution performed with DERIVE 6. 


4 2 2 2 

</(- z + 6 -z - 1 + 4 -l-z-(z - 1 )) + i • (— z - 1 ) 

W = - V 

2 + 2 • l -z 


4 2 2 2 

</(- z + 6-z - 1 + 4-L-Z-(Z - 1)) + L-(Z + 1) 

2 + 2 -l-z 


#5: SOLVE 


r Z — L W-L ^ 

: r- = °- z = 

V Z + L W + L ) 


4 2 2 2 

7(- w + 6-w - 1 + 4-l-w-Cw - 1)) + l • (— w - 1) 

z = - v 


2 + 2 • l - w 


4 2 2 2 

</(- w + 6-w - 1 + 4 -l-w-Cw - 1)) + l-(w + 1) 

2 + 2 -l-w 


Gert von Morze, Hannover, Germany 

As a Power User of DERIVE I came across this bug in Version 3.02: 

SUM ( COMB (n,k) , k , 0 , n ) = ?? 

Can you tell me why? 

DNL: As I have learned form the DERIVE manual DERIVE uses "antidifferences" to calculate 


sums: 


D 
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"As with antiderivatives , closed-form antidifferences may not exist in terms of the op- 
erators and functions known to DERIVE. Even when such an antidifference exists, there is 
no known method that is guaranteed to find it. " (DERIVE manual, page 202).. 

It could be a little consolation that in this case unlike to the problem above, even the Tl 92 
gives up. 

n n 

#6: I COMB(n p k) = 2 Solution performed with DERIVE 6. 

k=0 


Ian D 'Souza, Montreal, Canada (ian@vir.com) 

DERIVE Users: Help! I need answer for this one. 

Em running WINDOWS95 on a Pentium 120 with 32Mb RAM. DERIVE for WINDOWS can do the 
following integration without problems: 

=> Declare a , b as Real (positive) 

==> then Simplify 

INT(x A 4 EXP(-2a A 2*b A 2*x A 2 / (a A 2 + 3*b A 2)),x,0,inf) => works as expected. 

==> Now try to Simplify this one: 

INT(x A 4 EXP(-2a A 2 *b A 2 *x A 2 / (a A 2 + 3*b A 2)) / (a A 2 + 3*b A 2) A (7/2),x,0,inf) 

==> it gets stuck. 

==> change the exponent of x from 4 to 5 (no change to the power of x inside the exponential func- 
tion), Now it can do it! 

It gets stuck for all even powers of x > 4 and works for all odd powers of x. 

=> Substitute the denominator in the exponential function, (a A 2 + 3*Z? A 2), with G (declared Real - 
positive) and Simplify again - now it works no problem! ! Also seems to work when substitut- 
ing for the denominator of the expression itself. 

The previous revision DERIVE v4.00 gets stuck for all powers of x that are > 4 (even OR odd). 

My old DERIVE v2.01 can do all these types of integrals with ease. 

00 

( 7 2 

x n e x dx with k = positive are well defined. 

o 

*** What's going on? *** Why does the substitution with G work? 


dx = 


p!2 
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# 1 : 


# 2 : 


2 2 2 
4 - 2- a ■ 2- b -x 

x -e 


2 2 7/2 

(a + 3- b ) 

0 


2 2 2 
& - 2- a -2-b -x 

x -e 

dx 

2 2 7/2 

(a + 3- b ) 

0 


3-Jtt 


2 27/2 5 

256- (a + 3-b ) - | a- b| 


Solution performed with DERIVE 6. 


1 


662 27/2 

64- a ■ b ■ (a + 3- b ) 


Some of my other experiences with v4.00 point to a bug in the parser. Version 4.01 fixes some but 
others remain. DERIVE 4.00 and DERIVE XM have problems doing some integrals with very large 
(but constant) expressions of the form 

INT(BigConstant, 0, 0, 2n) 

where BigConstant is a very large expression, NOT a function of x. Where my old DERIVE v.2.01 
correctly returns 2*7i*Bigconstant, DERIVEXM and v4.00 seem to hang. I haven't fully experimented 
with v4.01, but it can now do some trigonometric integrals that v2.01 can do, but that XM and v4.00 
can't do. I've played with all the Manage Trig settings. 

One more thing: Note that v4.01 will not read a .MTH file (containing multicharacter Greek letters) 
created by v4.00! The extended ASCII characters used to represent concatenated Greek letters seems 
to have changed. 

Thanks in advance for any answers! Ian 

Llorens Fuster, Valencia, Spain (llorens@mat.upv.es) 

Dear DERIVERS, Let a function F(u,i):= 

(where u is for example a matrix, but that is not important). I want to programme the 'iterative - 
nested' function: 

F(u,l) 

F(F(u,1),2) 

F(F(F(u,1),2),3) 


F(F(F F(u,l),2),3) )n) 

where n (the number of iterations) depends on u. 
IS IT POSSIBLE? 


a* | 


I ^ I 
I a* I 


I ** i 
I a* I 


I ^ I 
I a* I 


I ^ i 
I a* I 


I ^ I 
I a* I 


I ^ i 
I a* \ 


I ** I 
I a* | 


I ^ \ 
I ^ \ 


I ^ I 
I | 


Another question: Is it possible to sort the elements of a vector? Thank you! 


DNL: I can only answer the second question: In DNL#1 3 you can find a sort routine: 

(The routine is on the next page. It is only of “historical” interest, because in the 
meanwhile the SORT-routine has been implemented, Josef) 
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The ERF(l) table (WIKI): 


0 . 771743 3 3 




0.82526296 

0.84310283 



0.83836778 

0.84273605 

0.84271160 


0.84161922 

0.84270304 

0.84270083 

0.84270066 

0.84243051 

0.84270093 

0.84270079 

0.84270079 



ERF(x) 



ERF(l) 

0.S4 2 7 0079 2949 714 
1 

2 2 

EXP (- t ) dt 

Jtt 

0 

0.842700792949714 


I found also a German Wiki-information, but according to my understanding there are some typos in it 
and following the instructions there the algorithm should not work (wrong subscripts?) 

I wrote a DERIVE program following this recursive procedure - which is based on the corrected pro- 
cedure provided on the German WIKI-page [3] . 



Then / n?n will deliver the n th approximation. 

The next diagram explains the calculation scheme: 


h,i 

h,i 



h,\ h,2 h,3 

U,\ ^ U,2 U,3 



Ua 


The program is displayed on the next page. 

Now I wanted to know - if my implementation is working and how the calculation times are perform- 
ing compared with DERIVE’s numerical integration, and with Terence’s and Steven’s procedures. 

I am starting with Terence’s integral from page 3. 

Let at first do DERIVE its job - needing nearly half of a second for performing the integration: 

[3] http://de.wikipedia.org/wiki/Romberg-Integration 
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4 

' SIN (x) 

#3 : dx 

. X 

0 

#4: 1. 7 5S 2 03138949 05 

0.469 sec 

This is my DERIVE Romberg program: 

rom_prog(f, x, a, b, n, eps, st, rtab, n_, k_, r_, i) := 

Prog 

st := VECTORCCb - a)/2 A (n_ - l)/2 • (LIM(f , x, a) + LIM(f, x, b) + 

2-I(LIM(f, x, a + i • (b - a)/2 A (n_ - 1)), i, 1, 2 A (n_ - 1) - 1)), n_, n) 

rtab := [APPEND([stil] , VECT0R(" ", j, n - 1))] 

n_ := 2 
Loop 

If n_ > n exit 
k_ := 2 

r_ := [st|n_] 

Loop 

If k_ > n_ exit 

i := rtatu(n_ - l)i(k_ - 1) + (rtatu(n_ - l)|(k_ - 1) - 
FI RST( REVERSE (r_)))/ (2 A (2 - 2-k_) - 1) 
r_ := APPEND(r_, [i]) 
k_ :+ 1 

rtab := APPEND(rtab, [APPEND(r_, VECT0RC" " , j, n - n_))]) 

If ABS(rtabin_in_ - rtabj,(n_ - l)i(n_ - 1)) < eps a 

ABS(rtabj,(n_ - 2)j,(n_ - 2) - rtabj,(n_ - l)i(n_ - 1)) < eps 
RETURN rtabj,n_ 4 ,n_ 
n_ :+ 1 
rtab 

Now let's compare: 

Terence Etchells: 0.078 sec 

5 8 12 

1.75804691765892 28.0139169547487 1792.45321974934 4.58861105626941-10 4.6987333062835-10 1.9 246 01049 3 0011 -10 

5 8 

1.75819500516204 28.0140294546536 1792.45365177833 4.58861112473565-10 4 . 69873331065501- 10 

5 

#6: 1.75820265343736 28.0140362330084 1792.45367852957 4.58861112900477-10 

1.75820310895234 28.0140366526776 1792.45368009761 

1.75820313707965 28.014 0366 788448 

. 1.7582031388323 

The first column contains the initial values which are found by Simpson’s approximation. According 
to Steven’s explanation, these values are the 2 nd step values in his procedure starting with approxima- 
tions obtained by applying the trapezoidal rule. Try finding the Simpson values in the next tables! 
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Steven Schonefeld: 0.047 sec 

' 1 4 1.62159875234603 *** *** *** 

2 2 1.72009680299869 1.75292948654958 *** *** 

4 1 1. 74855938899386 1.75804691765892 1.75838807973288 *** 

#5: S 0.5 1.75578610112 1.75819500516204 1.75820487766225 1.75820196969287 

16 0.25 1. 75759851535802 1.75820265343736 1.75820316332238 1.75820313611063 

32 0.125 1.75805196055376 1.75820310895234 1.75820313932001 1.75820313893902 

. 64 0.0625 1.75816534294818 1.75820313707965 1.75820313895481 1.75820313894901 

And this is my rom_prog: 

r sin go ^ 

#7: nom_pnog , x, 0, 4, 7 

' 1.62159875234603 
1.72009680299869 1.75292948654958 

1. 74S5593S8993S6 1.75804691765892 1. 75S3SS079732SS 

#8: 1.75578610112 1.75819500516204 1.75820487766225 1,75820196969287 

1. 75759851535802 1.75820265343736 1.75820316332238 1. 75820313611063 1. 7 58 2 0314 06848 2 

1. 75805196055376 1.75820310895234 1.75820313932001 1.75820313893902 1. 75820313895011 1. 75820313894841 
. 1,75816534294818 1.75820313707965 1.75820313895481 1, 75820313894901 1. 75820313894905 1.75820313894905 
My rom_prog: 0.047 sec 

f SIN(x) -8") 

#9: nom_pnog| , x, 0, 4, 7, 10 I = 1.75820313894905 

Including considering the required accuracy: 0.047 sec 

The last column of the matrix in #8 is not displayed. 

One can enter the required accuracy as last parameter and then the output is only the value of the inte- 
gral if the number of steps is sufficient enough. 

Finally I ask rom_prog to calculate the ERF(l)-value: 

(2 2 

#10: nom_pnog EXP(- x ), x, 0, 1, 6 

L Jn 

' 0.771743332258053 
0.825262955596749 0.843102830042981 

0.838367777441205 0.842736051389356 0.842711599479115 

# 11 : 

0.841619221244767 0.842703035845955 0.842700834809728 0.84270066394196 

0.842430505490232 0.842700933572054 0.84270079342046 0.842700792763488 0.84270079326867 

. 0.842633227681257 0.842700801744931 0.842700792956457 0.842700792949091 0.84270079294981 9 0.842700792949508 

0.016 sec 

f 22 -6^ 

#12: nom_pnog EXP(- x), x, 0, 1, B, 10 = 0.842700792949508 

L Jn ) 

0.046 sec 


This looks pretty good, doesn’t it? 
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POLIEDROS 1.0 


P.Familar Ramos 


POLIEDROS Version 1.0 is a PC 
program that represents graphically 
regular polyhedra (Platonic and 
Kepler-Poinsot solids) and semi- 
regular polyhedra and their duals 
(Archimedian and Catalan solids) 
through Aero Spin™ . The program 
includes a study of the metric prop- 
erties of each polyhedron (inform- 
ing about the number of faces, ver- 
tices and edges, the vertex configu- 
ration, the dihedral angles, the cir- 
cumradius, midradius and inradius, 
the surface area and volume). The calculations of these metric relationships have been avaluated with 
DERIVE ™. The user is allowed to select the program language: English or Spanish. 

Do you know a "Disdyakistriacontahedron"? And if so, do you also know that it is dual to a "Rhombi- 
truncated Icosidodecahedron"? And at last, do you how it looks like? POLIEDROS gives an answer: 


1. 

Platonic solids 

5. 

Dual relationships 

2. 

Kepler-Poinsot solids 

6 . 

Metric relationships 

3. 

Archimedean solids 

7. 

Language 

4 . 

Catalan solids 

Q. 

Quit 


a . 

Tr iak istetrahedron 

j- 

D isdyakisdodecahedron 

b. 

Tetrak ishexahedron 

k. 

Disdyakistriacontahedron 

c . 

Tr iak isoctahedron 

1 . 

Pentagonal icos itetrahedron 

d. 

Pentak isdodecahedron 

m . 

Pentagonal hexecontahedron 

e . 
f . 

Tr iak is icosahedron 
Rhombic dodecahedron 

n . 

Pentagonal dipyramid 

g- 

Rhombic tr iacontahedron 

o . 

Pentagonal deltohedron 

h. 

Deltoidal icos itetrahedron 
Deltoidal hexecontahedron 

Select option 






lUUUCLOIICUl UM 




Dua 1 

disdyakistriacontahedron 




Faces 

62 = 12 { 10 } *20 <6 > +30 




Uert ices 

120 




Edges 

180 




Uertex configuration 

(4.6.10) 




Dihedral angles 

{4}-{6> : 159" 5' 41" 





{4>-fl0}: 148° 16' 57" 





{6J--C10} : 142* 37’ 21" 




Circumradius 

3.80239 44998 51293 58476 




Midradius 

3.76937 71279 21716 60267 




Inrad i i 

3.73606 79774 99789 69640 

{4> 




3.66854 24806 72585 73611 

{6} 




3.44095 48011 77933 84551 

{10} 



Area 

174.29203 03423 23920 88293 



Volume 

206.80339 88749 89484 82045 



Press any key to continue... 









Very interesting for you is the fact, that POLIEDROS includes the ACROSPIN.EXE file. Math Ware 
has granted AWR Software a non-transferrable license agreement, subject to renegotiation, to use 
AcroSpin in POLIEDROS. 

Available at AWR Software, Huertos 21, 46500 SAGUNTO (Valencia), Spain, FAX 96 266 34 07 
(3700 Ptas) 

An Introduction to the Mathematics of Biology by Yeagers, Shonkwilerand Herod 

The authors of this textbook have adopted the philosophy that mathematical biology is not merely the 
intrusion of one science into another but has a unity of its own. The biology and mathematics are 
equal; they are complete and flow smoothly into and out of one another. 

The book has several important features that the authors have developed from their classroom ex- 
perience. A unique feature is the use of a CAS, Maple, in parts of every chapter. The models can eas- 
ily be transferred to other CA-systems as DERIVE. Graphic visualizations are provided for all the 
mathematical results. 

The chapters are: Biology, Mathematics, and a Mathematical Biology Laboratory; Some Mathematical Tools; 
Reproduction and the Drive for Survival; Interactions between Organisms and their Environment; Age- 
Dependent Population Structures; Random Movements in Space and Time; The Biological Disposition of Drugs 
and Inorganic Toxins; Neurophysiology; The Biochemistry of Cells; A Biomathematical Approach to HIV and 
AIDS; Genetics . 

The text has extensive exercises, problems, and examples, along with references for further study. 
Birkhauser Boston 1 996, ISBN 1 55953 1 07 X, 434 pages, Hardcover, DM 1 1 8.-, AS 862.-, sFr 98.- 
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This contribution from 1996 is of “historical and nostalgic” interest. In 1996 there 
was no 3D plot possible from within DERIVE. What we could do was producing 
3D-projections (isometric and other parallel projections, or perspective projec- 
tions applying self made projection procedures). And there was David Parker’s 
inexpensive ACROSPIN. We could call ACROSPIN from DERIVE - but only for 
functions of two variables - not for all other kinds of 3D objects. 


What I was demonstrating in 1996 can be reproduced in the DOS-environment 
even today - if you have ACROSPIN available 111 . But there is one obstacle. Com- 
pare the two files below. Both files have the same source and are saved as 
BASIC-files (what can be done in DERIVE 6 via the Write-option). 


— 


Load ► 


Write 

Basic File... 

II Handheld ► 

C File... 

Page Setup... 

Print Preview 

Print... Ctrl+P 

Fortran File... 

Pascal File... 

Rich T ext Format File . . . Ctrl+R 




This is how it looked like with DERIVE 3.14. 


#154:FIGUR_LKdian, diamw> 

1 

fflbb : 

i .7], [2.9, 0, i 

TRANSFER SAUE: Derive idAUd C Fortran Pascal Options State 
Enter option 

Approx<ttl54> D:\D0KUS\DNL\DNL24\ Free: 68* 

Derive Algebra 


Derive 3.14 


[[2,0.83,1.7], [2, -0.83, 1.7], [2. 9, 0,1. 2], [2,0.83,1.7], [2. 9, 0,1. 2], [3.7,- 
1.53,0], [4,0,0], [3.7,1.53,0], [2. 9, 0,1. 2], [4,0,0], [3. 7, -1.53,0], [2,0,- 
2. 2], [4, 0,0], [3. 7, 1.53,0], [2, 0,-2. 2], [0,0, -4]] 


Derive 6.10 


[ [2,0.83,1.7] , [2,-8. 3E-1 , 1.7] , [2. 9, 0,1. 2], [2,0.83,1.7], [2. 9, 0,1. 2], [3.7,- 
1.53,0], [4,0,0], [3.7,1.53,0], [2. 9, 0,1. 2], [4,0,0], [3. 7, -1.53,0], [2,0,- 
2.2] , [4,0,0] , [3.7,1.53,0] , [2, 0,-2. 2] , [0,0, -4] ] 


Compare -0.83 and -8.3E-1. I believe that it would be wasted time to change my BASIC- 
program ACD from 1996 and compile it again, because we do have now 3D-plots in 
DERIVE. 

In the following contribution I will add the respective DERIVE 6 plots. You may compare. 

Much fun reading or rereading DERIVE - ACD - ACROSPIN. 

Josef 

[11 l have a little present for all of you who have still the DERIVE DOS-version: David Parker 
gave permission to distribute the original ACROSPIN to the DUG-Members. you can 
download it among the other DNL24-files. Many thanks to David. 
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Interaction between DERIVE and ACROSPIN with ACD. 

Josef Bohm, Wtirmla, Austria 

Using DERIVE 3.x you can save surfaces as ACD-files if the surface is defined as a function of two 
variables f(x,y). For other 3D solids it is not possible to have a direct conversion to an ACD-file for 
running ACROSPIN-animations. It is also impossible to produce ACROSPIN-fUcs for space curves 
and for polyhedrons from the DERIVE-e nvironment. I must admit that I am a fan of ACROSPIN and I 
would like to recommend strongly to buy this piece of software. Nevertheless you may find this con- 
tribution also useful if you don't want to use ACROSPIN. Using the functions ISOMETRIC and 
COPROJECTION from DERIVE' s utility file GRAPHICS.MTH you can also produce impressive 
plots. See the end of this paper. Fortunately there are a lot of other DERIVIANS which like geometry, 
so I am glad to announce for the next DNLs some contributions to produce several mappings of ob- 
jects (H. Kummel a.o.) and - a hidden line algorithm for DERIVE from our friend Hubert Weller. 

It is easy analysing the format of an ACD-file and then to edit a few 3D-points for simple objects for 
one's own ACD-file. To do so for objects consisting of many points this will be a boring work. 

So I remembered my "programming past” and wrote a tool, ACD.EXE which can help. I am sure that 
DERIVE and ACROSPIN as well may benefit 

Now you are able to animate space curves, polyhedrons and surfaces (given in any parameter form) in 
any combination of objects, colours and layers. Using an algorithm of Richard Schorn from Kauf- 
beuren you can produce analglyphs (red-green-pictures) to obtain a stereographic presentation of the 
object. At this place I want to thank Mr Schorn for his comments and support. Many letters were ex- 
changed between Kaufbeuren and Wiirmla to share experiences. 

First you have to produce the list of 3D-points of the object in DERIVE. Work in Approximate Mode 
and use 3 digits in Notation (from the Options submenu). 

If you connect the points then they will form a space curve or a polyhedron or only a polygon in 
space. Save the DERIVE-expression as a BASIC-file. So the most important thing is to create lists of 
points in the right order. These lists of points can be used for many other projections (parallel or cen- 
tral perspective, ) or the projection process can be done by another tool - like ACROSPIN. You 

can also use this lists in connection with ISOMETRIC. 

So I will start with a toolbox for producing this lists of points for 3D objects. 

"P3D . MTH" 

[ InputMode : =Word, Precision : ^Approximate ] 

[ PrecisionDigits : = 4 , Notation : ^Decimal , NotationDigits : =3 ] 

"The next 7 functions are from:" 

"File GRAPHICS .MTH, copyright (c) 1990 by Soft Warehouse, Inc." 

COPROJECTION (v) : =VECTOR (VECTOR (u_ SUB n_, u_, v) , n_, DIMENSION (v SUB 1)) 

SPHERE (r, theta, phi) :=r* [SIN (phi) *COS (theta) , SIN (phi) *SIN (theta) , COS (phi) ] 
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CYLINDER (r, theta, z) := [r*COS (theta) , r*SIN (theta) , z] 

CONE (alpha, theta, z) := [z*SIN (alpha) *COS (theta) , z*SIN (alpha) *SIN (theta) , z] 
NORMAL_VECTOR(v, t) :=SIGN (DIF (v, t, 2) ) 

BINORMAL (v, t) :=SIGN (CROSS (DIF (v, t) , DIF (v, t, 2) ) ) 

SPACE_TUBE (v, t, r,phi) :=v+r* (SIN (phi) *NORMAL_VECTOR (v, t) +COS (phi) * 

BINORMAL (v, t) ) 

"The following is not part of GRAPHICS . MTH : " 

FIG (v_, pl,p2,pla,ple,n,p2a,p2e,m) : =VECTOR (VECTOR (v_, p2,p2a,p2e, 

(p2e-p2a) /m) ,pl,pla,ple, (ple-pla) /n) 

RO(obj,n) :=VECTOR(obj . ROTATE_Z (2*pi*k/n) , k, 1, n) 

FIGURJaT (p_, w ) : =APPEND (VECTOR ( [ELEMENT (p_, ELEMENT (w_, kl_) ) ] , kl_, 

DIMENSION ( w ) ) ) 

FIGUR_E (p_,e_) :=APPEND (VECTOR ( [ELEMENT (p_, ELEMENT (e_, kl_, 1 ) ) , ELEMENT (p_, 
ELEMENT (e_, kl_, 2) ) , [inf, inf, inf] ] , kl_, DIMENSION (e_) ) ) 

PI (m, zn) :=VECTOR (VECTOR ( [ELEMENT (m, i, 1) *COS (phi_) , ELEMENT (m, i, 1) *SIN (phi_) , 
ELEMENT (m, i, 2) ] , phi_, 0 , 2 *pi , 2 *pi/zn) , i , 1 , DIMENSION (m) ) 

P2 (m, zn) : =COPROJECTION (PI (m, zn) ) 

REV (m, zn) :=APPEND(P1 (m, zn) , P2 (m, zn) ) 

You can create a polyhedron with a list containing all the edges of the solid: 

[[P1,P2],[P2,P3], [Pi,Pj], ], with Pi = [xi,yi,zi]. Save this file as a BASIC-file. I prepared some 

DERIVE - functions to produce polyhedrons from a list of its points and either a list describing the 
way how to connect the points or a list (matrix) containing the edges: P_WAY(points,way) and 
P_EDG(points, edges). 

points : = [ [0,0,0] , [3, -2, -1] , [1, 5, 10] , [1, 2, 5] 
way : = [ 1 , 2, 3, 1,4, 3, 4, 2] 

edges : = [ [1,2] , [2, 3] , [3, 1] , [1, 4] , [2, 4] , [3, 4] ] 

P_WAY (points, way) 

P_EDG (points, edges) 

Both functions will when Simplified return the points of a pyramid. 

Or you have produced a family of parameter lines with parameters u and v using the powerful 

VECTOR(VECTOR(f(u,v), ^-construction. Using COPROJECTION from GRAPHIC S . MTH you 

will have both families of parameter lines with either worv constant. Save the two expressions as dif- 
ferent BASIC-files, leave DERIVE and then call ACD. The desire to produce ACROSPIN- 
demonstrations of solids of revolutions was the inspiring idea to create a tool like ACD to make that 
possible. 

Executing ACD, you will find a menu in which you are asked to enter the type of object which you 
want to prepare for ACROSPIN. You are able to combine several different objects in one ACD-file 
with both individual colours and layers. See now some examples: 
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From DERIUE to ACROSPIN 

Which kind of object do you want to display? Make your choice: 


Space Curve or Polyhedron (from a list of points) ... 1 
in Stereo - Uision 4 

Surface (shou the families of parameter lines) 2 

in Stereo - Uision 5 

or shou the complete surface 3 

in Stereo - Uision 6 

Polyhedron (from a list of edges) 2 

in Stereo - Uision 5 

Discrete points -1 

in Stereo - Uision 7 

Quit 0 


Your choice please: 


From my point of view the most interesting thing is the fact that I can interact between the wonderful 
capabilities of DERIVE and the demonstrating facilities of ACROSPIN. We could use this tool to train 
the students' imagination of 3D space. Give them the task to design 3D objects by the coordinates of 
the vertices and the edges. Let them construct well known geometric solids like tetrahedrons, octahe- 
drons, pyramids, parts of a cube, and so on. They immediately can see if they are right or not. 

Rotating in ACROSPIN you will obtain top- front- and side view of your object. You have to imagine 
all the following pictures in different colours and layers. So you can switch off and on different parts 
of the objects. See a composition of the objects axes, house and tower. 
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In DERIVE 6 you have to Insert (F4) the objects 
separately 

(FIGUR_W(axes , axesw), 

FIGUR_W(tow, toww), 

and FIGUR_W(house , housew) respectively. 

Then change the Scheme in the Color Plot Window 
to Custom and choose the color of your choice for 
the Grid. 




I called the first combination of objects "VILLAGE". In the second figure you can see the eighth of a 
diamond. Rotating this part seven times using DERIVE the whole figure will emerge: 

(You can find the objects in P 3D.MTH). 

"The Diamond - (the girls' best friend)" 


diam: = [ [2, -0.83, 1.7], [2. 9, 0,1. 2], [3.7, -1.53,0], [4,0,0], [2, 0,-2. 2], 
[0, 0, -4] , [2,0.83,1.7] , [3.7,1.53,0] ] 

diamw : = [7, 1,2, 7, 2, 3, 4, 8, 2, 4, 3, 5, 4, 8, 5, 6] 

FIGUR_W (diam, diamw) 

RO(FIGUR_W (diam, diamw) ,8) 


Save both results as different BAS-files, then run ACD. For the first part apply option 1 and then ap- 
pend the whole object in another colour using option 2. 




RO(FlGUR_W(di am, diamw) ,8) results in the left plot in DERIVE 6, you have to 
APPEND(RO(FIGUR_W(di am, diamw), 8)) 
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The left picture is an analglyph of a C60-molecule model (see the contribution submitted by Richard 
Schorn in DNL#21). You have to imagine red and green lines and viewing them through a red-green- 
glass. 



The right picture is composed from a sphere and a space curve - a loxodrome. You can make visible 
the single points applying option 1 and then produce the closed curve using option 2. I’ll lead you the 
way how to do it. 



At first we have to produce the list of points using DERIVE : 

#2 8: [ c : =0 . 2 , rl : =#e A (c*phi) ] 

# 29 : [ rO : =1/ SQRT (l+rl*rl) ,w:=ATAN(rl) , r : =rO*SIN (w) ] 

#30 : VECTOR ( [r*COS(phi) ,r*SIN(phi) ,rO*COS(w) ] ,phi, -20, 0, 0 . 25) 

#31 : VECTOR ( [r*COS(phi) ,r*SIN(phi) ,rO*COS(w) ] , phi, 0, 20, 0 . 25) 

#32 : VECTOR (VECTOR ( [0 . 5*SIN (phi) *COS (theta) , 0 . 5*SIN (phi) *SIN (theta) , 

0 . 5* COS (phi) +0.5] , phi, 0,2*pi,pi/10) , theta, 0,2*pi,pi/10) 

#33 : FIG ( [0 . 5*SIN (phi) *COS (theta) , 0 . 5*SIN (phi) *SIN (theta) , 0 . 5*COS (phi) +0.5], 
theta, phi, 0,2*pi,20,0,2*pi,20) 

#30 and #31 give the points of the loxodrome. I had to divide into two parts, because ACD can work 
only with vectors consisting of 100 elements (each of which can contain again 100 components). 
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Expression #33 is the parameter form of a sphere. Using my FIG-function you can avoid the bulky 
VECTOR (VECTOR ( ) ) - command. 

Approximate #30, #31 and #33. 

Save each of the results as a single BAS-file with different names, eg LOX1, LOX2 and SPH respec- 
tively. Quit DERIVE. Take care that ACD.EXE and ACROSPIN.EXE are in the same directory. Then 
call ACD. 

Choose option 1 , because we want to see the points which build up the space curve, 


LOX1 

[ENTER] 

the name of the first BAS-file 

LOXO [ ENTER] 

that's the name of the ACD-file to be created now 

15 

[ENTER] 

the points should be white (colour code 15) 

1 

[ENTER] 

the 1st layer 

y 


because we have not finished, repeat with LOXl, you will not be 
asked for the ACD-file 's name once more. If you would finish now you 
would obtain only the points (next figure!!), but we will continue: 

l 

[ENTER] 

we will see the space curve 

LOXl 

[ENTER] 

the same point list will be used 

4 

[ENTER] 

we want to have a red curve 

2 

[ENTER] 

it should be another layer 

y 

[ENTER] 

because we are not yet ready, repeat the two steps with LOX2. 

If we would finish now we 

could see the points and the curve in two layers, which could be toggled on 

and off using 

the [ ENTER] 

-key in combination with the [ 1] - and[ 2] -key. As we want to add the 

sphere we go on once more 

and proceed: 

3 

[ENTER] 

the whole surface of the sphere 

SPH 

[ENTER] 

the name of the corresponding BAS-file 

14 

[ENTER] 

let 's have a yellow sphere 

3 

[ENTER] 

in the 3rd layer 

n 


we have finished. Now start typing ACROSPIN LOXO [ ENTER] 
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These are the DERIVE 6 commands: 


APPEHD(C0PR0DECTI0N(FIG([0. 5- SIN(tfi) -C0S(6) , □. 5- SIN 0*0 ■ SIN( 6 ) , 0.5.C0S(<|>) +0.5], 6 , $, 0, 2-n, 20, 0, 2-n, 20])) 
APPEND(FIG([0. 5- SIN(t|>) -005(9) , 0. 5- SIN(<|>) ■ SIN(B) , 0.5.C05((|)) +0.5], 0, 4 ., 0, 2-n, 20, 0, 2-n, 20)) 

VECTOR ([r-COS(tfi) , r-SINOtO, r0.C0S(»], 4 ., -20, 20, 0.25) 

VECT0R([[r-C0S(4i) , r.SIN(t|j), rO-COS(w)]], (|), -20, 20, 0.25) 



Inspired by books dealing with other CAS - packages I tried to produce a more sophisticated anima- 
tion: an elliptic torus with a torus knot line on it. And to make the space curve more impressive I su- 
perimposed its tube. (Using the utility file GRAPHICS.MTH from SWHH). 


#103 : TKN (a, b, c, p, q, t) := [ (a + b • COS (q • t ) ) • COS (p • t ) , 

(a + b-COS (q-t) ) -SIN(p-t) , c-SIN(q-t)] 

#104 : TKN (8, 3, 5, 2, 5, t) 

#105: [3-COS(2-t) -COS(5-t) + 8-COS(2-t), 3 • SIN (2 - 1 ) - COS ( 5 - 1 ) + 8-SIN(2-t), 
5 • SIN (5 • t) ] 


n -1 


#106 : VECTOR 


TKN (8, 3, 5, 


#107 : ELL TOR (a, b, c, 0, 


#108 : FIG (ELL_TOR (8, 3, 5, 

#109 : SPACE_TUBE (TKN (8, 3, 

# 1 1 0 : FIG ( SPACE_TUBE ( TKN ( 8 
2-n, 20) 


2, 5, t) , t, 0, 2-n, 

40 - 1 

<t>) := [ (a + b-COS (<D) ) -COS (0) , a + b • COS (<D) • SIN (0) 

C • SIN (4>) ] 

0, <t>) , 0, <t>, 0, 2-n, 40, 0, 2-n, 40) 

5, 2, 5, t) , t, 2, 4>) 

, 3, 5, 2, 5, t) , t, 1, <t>) , t, <t>, 0, 2-n, 80, 0, 
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You have to approximate the expressions #106, 
#108 and #110, save them in three different 
BAS-files and then use ACD. Please be patient 
approximating expression #110. Using ACD I 
chose colour 1 for the space curve, colour 4 for 
the torus and colour 14 for the tube. 


DERIVE 6 graphs are below. 
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If you know the parameter form of a surface it is very easy to produce the family of parameter lines 

with a VECTOR(VECTOR( )) function and then to apply option 2 or 3 from ACD.EXE to create 

the animation of the surface. You all will know the famous Moebius Strip. Its parameter form is given 
by: 


F(u, v) 




With a = 1 we obtain the DERIVE representation. 

VECTOR (VECTOR ( [COS (u) +v*COS (u/2 ) *COS (u) , SIN (u) +v*COS (u/2 ) *SIN(u) , 
v*SIN(u/2) ] ,u,0,2*pi,pi/20) ,v, -0.3, 0.3, 0.1) 

instead of VECTOR (VECTOR ( [ ) ) you can use my FIG ( ) -function: 

FIG ( [COS (u) +v*COS (u/2) *COS (u) , SIN (u) +v*COS (u/2 ) *SIN (u) , v*SIN (u/2 ) ] , u, v, 
0, 2 ^pi, 40, -0.3, 0.3, 60) 



#111 : FIG (mb, u, v, 0, 2*n, 40, -0.3, 0.3, 6) 

r d d ■ 

#112: (lim mb) + t- CROSS lim — mb, lim — mb 
v->0 L v ->0 d u v->0 dv 


#113: [t-COS (u) -SIN(0.5-u) + COS (u) , t • SIN (u) • SIN ( 0 . 5 • u) + SIN(u), -t • COS ( 0 . 5 • u) ] 
#114 : M v = 0 gives the circle in the middle of the strip" 

#115: [0-COS (u) -SIN(0.5-u) + COS (u) , 0 • SIN (u) • SIN ( 0 . 5 • u) + SIN(u), -0 • COS ( 0 . 5 • u) ] 
#116 :"t = 0 --> pedal points of the normals" 

#117 : [COS (u) , SIN (u) , 0] 

#118 : "t = 0.5 --> end points of the normals" 

#119: [0.5-COS (u) -SIN(0.5-u) + COS (u) , 0 . 5 • SIN (u) • SIN ( 0 . 5 • u) + SIN(u), 

-5-COS (0.5-u) ] 

rr COS (u) SIN (u) 

#120 : VECTOR 

-L 0.5-COS (u) -SIN (0.5-u) + COS (u) 0 . 5 • SIN (u) • SIN ( 0 . 5 • u) + SIN(u) 


0 



0.5-COS (0.5-u) 


, u, 0, 2 • n. 


20 
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As you can imagine you can use this interaction between DERIVE and ACROSPIN not only to pro- 
duce ’’nice” pictures. In a comfortable way you can make visible results from differential geometry. So 
you see on the other Moebius strip the normals. You could add the tangents, the normals and binor- 
mals to space curves, the tangent planes and, and, 

Approximate expressions #111 and #120 and save them as you have done before. Then apply option 3 
for the strip and option 2 for family of normals (try option 3, then you will obtain a second strip! !??). 


There is another idea to use this interaction: make 
clear the meaning of the parameters. Let me dem- 
onstrate my idea: 

We create a spindle torus and then try to produce 
certain sections: 




Which values for the parameters are responsible for the various torus parts? 

Unfortunately I cannot print in colours - maybe in some years we will have a colour print DNL??? 

(Yes, now we have!!) 

As I have mentioned before, I wanted to demonstrate solids of revolution. I also wanted to make visi- 
ble the two families of parameter lines. So I produced a DERIVE function to create these surfaces from 
any given profile. We can design a profile by a list of points or by a - piecewise defined - function. I 
show one example for each possibility (left and right plot windows): 
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The next plots show the solids: you 
can see the profile together with the 
parallel circles, then the longitudinal 
intersections and at last the full grid. 
It is very nice to overlay the different 
layers and then to rotate the glasses 
round their axes, to zoom in and out, 
to translate, to accelerate, 
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In this collection you will find a rotated lemniscate, three quarters of a pseudosphere and another ro- 
tated curve. There exists a simple formula for producing a solid of revolution using any curve given in 
parameter form [x(t), y(t)]: 


F(t,<p) = [x(t) cos (p, x(t) sin <p, y(t)] 












The next example uses the zooming abilities of ACROSPIN: 

I produced the Henneberg Surface in different scales: 

I took another scale and zoomed in. Observe the dark spot in 
the centre of the graph at the right and what is hidden in it: 


2 

2 ■ SINH(s) ■ COS(t) SINH(3 • s) -COS (3 ■ t) 7 

3 

2 

2 * SINH(s) * SIN(t) + SINH(3 *s) *SIN(3 • t) , 

3 



2 ♦C0SH(2*s) « CGS(2 - 1) 
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It is obvious that it cannot be too difficult to produce presentations 
of intersecting surfaces together with the intersection curves, tan- 
gents, normals, binormals, osculating planes, You can also add 

labels as you will see in one of the next 
examples. 

I teach in a secondary school and I have some models of the conics 
sections. But using a tool like DERIVE it is convenient to produce a 
double cone, the intersecting plane, then to calculate the intersecting 
curve as a space curve and a bit more ambitious to shade the inter- 
secting surface. All this is hard to do by hands only, but let DERIVE 
do the calculations. The students have to know the strategy to obtain 
results which can be used to be represented by ACROSPIN. 


You can turn the model round and observe it from all directions. The different colours enforce the 
imagination. Switch off and on the layers. Compare the hyperbolic section with the parabolic and the 
other ones ... 

I want to represent a double cone, the intersecting planes, the inter- . 

section curves and I will try to add a shading for the intersection 2H=4 \ 

figure to make the picture more impressive. \ 

You can find the whole calculation in CONICS.MTH. The various ph 

CONIC*. ACD files show the different conic sections. 


I want to show the start of the ’’parabolic” part to give some com- 
ments on it. 


#3: cone := [3 • t • COS ($) , 3 • t ■ SIN ($) , 4 - 4-t] 


- 


n - 

" 

VECTOR cone, 0, 0, 

2 * n. 

r t , 

-1, 1, 2 

- 


20 - 1 

- 1 

[-1, 0, 4] + u- [3, 

l 

o 

+ v- [0, 

1 , 0] 



#54: sc 3 := 


V([3, 0, -4] • [ 3 , 0, -4]) 


r r 

sc 3 - 


1 

VECTOR pp, u, 0, 1, 

2 - 

, v. 

-3, 3, 0.5 


#57: SOLVE (pp = cone, [t, u, v] ) 

r 1 1 

#58: t = u = 

L 3 - (1 - COS (4>) ) 3 - (1 - COS (4>) 

#60: [ 3 • t • COS (<£) , 3 • t • SIN (<£) , 4 - 4-t] 


SIN (<£) 


1 - COS (4>) 


1 1 

#61: 3 COS (4>) , 3 SIN(<£), 4 

3- (1 - COS(<f)) 3 • ( 1 - COS (4>) ) 


3 • (1 - COS (4>) ) J 

COS (4>) SIN (4>) 4 • (3 -COS (4>) - 2) 

f r 

1 - COS (4>) 1 - COS (4>) 3 • (COS (4>) - 1) 
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#64: SOLVE 


4- (3-COS (O) - 2) 


= 0,0 


3- (COS (O) - 1) 


#65: 


O = ASIN 


-, O = ASIN 


3 • n n 

+ , O = - ASIN 

2 2 


r 2 


L 3 JJ 


#66: [O = -0.841, O = 5.44, O = 0.841] 

COS(O) SIN (O) 4 • ( 3 • COS (O) - 2) 


#67: VECTOR 


- L 1 - COS(O) 1 - COS(O) 3 • (COS (O) - 1) 

5.44 - 0.841 


, O, 0.841, 5.44, 


40 


#84: VECTOR 


- V (2 -w + 1) • SIGN (w + 1) 


V (2 -w + 1) • SIGN (w + 1) 


4 • (2 - w) 


4 • (2 - w) 


, w, -0.5, 2, 0.05 


In line #3 cone is the parameter form of the double cone, it is intersected by the plane pp in #52. 
Approximating #4 returns the family of lines building the cone, approximating #55 gives the grid of 
the plane. #54 is necessary to obtain a scaling factor to have a 1 by 1 grid. We intersect cone and plane 
in line #57, substitute for t in the cone’s parameter form and find a "space curve” in #63, which is our 
conic section - a parabola. As we don’t want to plot a parabola coming from anywhere and leaving for 
anywhere we have to find the parameter values which represent the intersection points of the parabola 
with the base. 

The shading should run normal to the parabola’s axis, so we reparametrize the curve to end with the 
shading’s vector in #84. We then approximate #4, #55, #67 and #84, save in different BAS-files, run 
ACD, and so on. 

The next plot shows one of the favourite examples of Bert Waits. So it is an honour for me to dedicate 
this DERIVE - ACD - ACROSPIN - product Bert and his "Power of Visualization". It is a real repre- 
sentation of the complex roots of a 5th order equation. I used DERIVE to find the modulo surface of 

z 5 6z 3 + 25z 0. 

The peaks are in the complex plane at the positions of the five solutions. I added the complex plane, 
the axes and their names "R" and "I" - both letters as special "space curves", given by a list of points. 
(Look at the file BERT.MTH on the diskette!!). The "R" can be seen in this position. The other graph 
shows another artificial range: a "fractal landscape" generated by a recursive algorithm in DERIVE , 
which I will present in one of the next DNLs. 

Bert's Complex Range 


The Fractal Range 
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|i ■ i ii 

From fractals it is not very far to chaotic behaviour. Josef Lechner was the first - besides Richard 
Schom - who checked ACD. He immediately tried successfully to visualize and to animate dynamic 
systems of three variables. During a phone call he mentioned the possibility now to animate the 
LORENZ attractor. So I will finish my interaction between DERIVE and ACROSPIN showing the 
famous LORENZ-attractor together with a zoom in from another "look out”. The list of points was 
generated with Josef Lechner’ s DERIVE - tool INTEGRAPH. (Proceedings of the DERIVE Days 
Diisseldorf). 

The following system of differential equations is the base of this attractor: 



At last I want to give an impression how to work without ACROSPIN . I produce a “static” representa- 
tion of the parabolic intersection of the double cone. As you can learn from the screen shot the trick is 
to to include ISOMETRIC and COPROJECTION from GRAPHICS.MTH at the appropriate place. 
Please compare with the according file on pages 31 and 32. Don’t forget to use COPROJECTION, 
other wise you would see only one family of parameter lines. 


? : 

J 

#85: UECTOR UECTOR ISOMETR IC (cone ) j 5, 

x^v\V\ 1 \TTJ77x 

L L 

C\\ \ \ \\ \l III/// / /\ 


(\\\\ \ \\ \\ ill//// / x\ 


VN\\\ \ \ \ \ \ 1 \I//Z// Z/XJ 

r r r 

mww\i I \ i j n / / ///xz 

#97 : COPROJECT ION UECTOR UECTOR ISOMETR 


#90: UECTOR juECTOR ^ISOMETR IC (eb3 ) , u, 0 


#98 : COPROJECTION juECTOR juECTOR [iSOMETR 

BPEllHH 

r r C0SC5) 


#93: UECTOR ISOMETRIC , — 

!* ' lx : * * S 

L L 1 - COSCS) 1 

1 

^^JuECTOR ISOMETRIC u, - J(Z u + 1) 







p34 


Josef Bohm: DERIVE - ACD - ACROSPIN 


D-N-L#24 


Among the downloadable files you can find in <ACD> all the MTH-files mentioned in this contribu- 
tion accompanied by ACD.EXE and a self extracting compressed file ACDZIP.EXE containing a lot 
of ready made ACD-files. And you will also find ACROSPINW 

I would like to ask you to produce your own ACDs. And to enforce this "CALL for ACDs" I 
invite you for a competition 


The ACD of the Year. 

The "ACD of the Year" will win one year free DUG membership. Deadline is 31 May 1997. Much 
luck. 

I am glad to announce another highlight for the next DNL. Our friend Sergey Biryukow from Moscow 
has produced a DERIVE tool to produce plots of implicit 3D functions. He saw ACD at Bonn, so he 
took care that his output is ACD compatible. 

Last question: Which object is hidden in this star? It is the Bottle of KleinW 
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Abstract 

In this work the general solution, or a particular one, to linear ordinary differential equa- 
tions, homogeneous or non-homogeneous, of third order with constant coefficients, is 
given, using DERIVE. A file with the appropriate DERIVE functions is included. 

1. INTRODUCTION 

DERIVE cannot be considered a programming language. However, new functions that make use 
of the operations and functions included in the program can be defined. In the program handbook it is 
said that DERIVE is not more than a collection of mutually recursive functions. 

We believe that (with the available information) differential equations have not been explored 
adequately. 

In the hand book, we are informed that DERIVE contains three files ODE1.MTH, ODE2.MTH 
and ODEAPPR.MTH to solve respectively ordinary differential equations of first, second-order and 
numerically. We have designed a file LODE3.MTH, to find the general solution, or a particular one, of 
an homogeneous or non-homogeneous linear differential equation of third-order with constant coeffi- 
cients: 

y'”+ ay” + by' + cy = f(x) 

2. DESCRIPTION OF THE FILE 

Firstly the roots of the characteristic polynomial x 2 + ax 2 +bx + c are determined. Considering 
that these roots appear in vectorial form where each element is an equation in the form x = “root”, the 
RHS function that selects the right member of the equation has been used. This function RHS is only 
present in DERIVE 2.58 and later. 

Next, the function F(a,b,c) has been defined to control the different kinds of roots that the charac- 
teristic polynomial can have. Thus, if F (a,b,c) > 0, the polynomial has a real root and two conjugate 
complex ones. If V(a,b,c) < 0, the polynomial presents three different real roots; and if F (a,b,c) = 0 it 
has only one real root with multiplicity three if and only if 3 b = a 2 . 

When the polynomial presents two real roots, one simple and the other of double multiplicity, a 
function P(a,b,c) has been defined, due to the fact that the SOLVE-function gives a vector with two 
components (one of them being the double multiplicity root). The sign of P (a,b,c) together with the 
relative magnitude of the vector components determine the position of the multiple/simple root in the 
vector. Taking into account the previous considerations about the roots, the function 
SGH(a,b,c,x,cl,c2,c3) calculates the general solution of the homogeneous equation. 
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As Y1 (a,b,c,x), Y2 (a,b,c,x) and Y3 (a,b,c,x) functions, defined the program, are linear independ- 
ent solutions of the homogeneous linear differential equation, it can be assured that if the parameter 
variation methods for finding a particular solution of the non-homogeneous equation has been used, 
then the corresponding equation system has only one solution which can be obtained by means of 
Cramer's rule. 


Then a function SGC(a,b,c,x,f,cl,c2,c3) has been defined. It calculates the general solution of 
non-homogeneous differential equations. 

Finally, the functions SPH(a,b,c,x,x0,y01,y02,y03) and SPC(a,b,c,x,f,x0,ly01,y02,y03) have 
been defined to calculate a particular solution to both homogeneous and non-homogeneous equations, 
respectively. 

It is just necessary to remember the names and arguments of the functions SGH, SGC, SPH and 
SPC in their right order. 


3. FILE LISTING 

"File LODE3.MTH" 

RPC (a,b, c) : =SOLVE (x A 3+a*x A 2+b*x+c, x) 

RCP (a, b, c) : =VECTOR (RHS (ELEMENT (RPC (a, b, c) , k) ) , k, 1,3) 
T1 (a,b, c,x) : =EXP (ELEMENT (RCP (a,b f c) , 1) *x) 

T2 (a,b, c, x) : =EXP (ELEMENT (RCP (a, b, c) , 2) *x) 


(As the file takes more than two pages and fortunately I can add the Diskette of the Year to 
this issue you can find the file LODE3.MTH on the diskette in subdirectory <MTH24>. So we 
will go on immediately to the examples. The paper submitted had included 8 pages with 81 
examples. I'll try to give a selection of some typical examples. Josef) 

4 . EXAMPLES 

A 486 DX4/75 has been used in the solutions of the differential equations below. The file has 
been optimized in order to get 

a) the solution to all third-order linear differential equations with constants of [1],[2] and 

b) its minimum execution time. 

I changed the original file for its use with later DERIVE versions. The first two expressions 
had to be adapted because of distinguishing between SOLVE and SOLUTIONS. 

Josef (2010). The first expressions read now: 

#1: File L0DE3_V6.MTH 

3 2 

#2: RCP(a, b, c) := S0LUTIONS(x + ax + b-x + c, x) 

Tl(a, b, c, x) := EXP((RCP(a, b, c)) -x) 

#3: 1 


# 4 : 


T2(a, b, c, x) := EXP((RCP(a, b, c)) -x) 

2 
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1) y" r - 2 y" - 5 y' +6y = 0; Check the solution! 


3-x x - 2-x 

#33: 5GH(-2, -B , 6) = c3-e + cl-e + c2-e 

3-x x - 2-x 

#34: y(x) := c3-e + cl-e + c2-e 

#3B: y rrr (x) - 2-y rr (x) - B-y r (x) + 6-y(x) = 0 


2) y m - y" = - 2y ; Check the solution! 

x -x 

#37: SGHC-1, 0, 2) = e -(c2-CGS(x) + c3-5IN(x)) + cl-e 

x -x 

#3S : y(x) := e -(c2-C05(x) + c3-5IN(x)) + cl-e 

#39: y rrr (x) - y rr (x) + 2-y(x) = 0 

,,f + 5 yff + 

3 ) y = — 

11 

x — 3 ■ x 

#41: SGH(5, 5, -11) = cl-e +e ■ (c2.C0S(y2-x) + c3- SINGfc.x)) 


4) y m - 4 y" + 7 y' = 6 y; ><0) = 1, y'(0) = y"(0) = 0 

Find the solution and check the result! 

2-x x 

#43: SPHC-4 , 7 , -6, x , 0, 1, 0, 0) = e - -/2-e -SIN(72-x) 

2-x x 

#44: y (x) := e - /2-e -SINO/2-x) 

#4B : y rrr (x) - 4-y rr (x) + 7-y r (x) - 6-y(x) = 0 

#46: [y (0) , y r (0), y rr (0)] = [1, 0, 0] 


5) y m - 3y" + 4y = x e 2x - cosx; Check the result! 


2-x 

#4S : SGC(-3, 0, 4, x, x-e - C0S(x)) 


2-x 3 2 

e -(9-x - 9-x + 6-x- (27-c2 + 1) + 162-cl - 2) -x 7-C0S(x) SIN(x) 

#49: + c3-e + 


162 


B0 


B0 


2-x 3 2 

e -(9-x -9-x + 6-x- (27-c2 + 1) + 162-cl - 2) -x 7-C0S(x) SIN(x) 

#B0: y(x) := ■+■ c3-e - + 


162 


B0 


B0 


#S1: y rrr (x) - 3-y rr (x) + 4-y(x) 


2-x 2 2-x 2 

e -(9-x - lB-x + 2- (27-c2 - 2)) e -(9-x + 12-x + 2- (27- c2 - 2)) 

#B2 : + C0S(x) 


27 


27 


#B3: needs resi'mpli'fyi'ng #B2 


■ 


#S4 : 
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6) y m - 2 y" + 5 y' = - 24e 3x ; y(0) = 4, /(O) = -1, y"(0) = 5 ; Check the initial conditions! 

3-x 


#56: 

spec- 

-2, 5, 0, x, 

- 24 -e 

, 0, 

4, -1, 

5) 




#57: 

SPC(- 

-2, 5, 0, x, 

3-x 

- 24 -e 

, 0, 

4, -1, 

5) = - 

3-x 
- e 

X 

+ e ■ (2-SINC2-X 

) - 2-COS(2-x)) + 7 

#58 : 

y(x) 

3-x 

:= - e + 

X 

e ■ (2- SIN( 

2-x) 

- 2 -CO 

S(2-x); 

) + 7 



#59: 

[y(o: 

), y’CO), y' 

'(0)] = [4 

, -1 

, 5] 






y w - 

2/' - 3/ 

+ 10 J = 


(34x - 

16) - 

10x 2 

+ 6x +34 


1 ) 

y (0) = 

3,y'(0) = y 

" (0) = 0 







#61: 

: SPC(- 

-2, -3, 10, 

- 2-x 

x | e 

(34- 

X - 16) 

- 10- : 

2 

x +6 

■x + 34, 0, 3, C 

), 0) 


~ £. ■ X £ 

#62 : x -e - x + 3 


8) Find the general solutions: 


/' + 5/ + 26/ - 150;; = 20e~ x 
y"’ + 5 y" + 26 y - 150j = 600 x 3 
4 


m At 

y + 4y 


tan 2x 


#64: 

SGC(5 , 

-x 

26, -ISO, x, 20-a ) 


#65: 

3- 

cl-e 

-x 

x 5-e - 4-x 

+ e -(c2-C0S(^: 

43 

34 -x) + c3-SIN(/34-x)) 

#66: 

SGC(5, 

3 

26, -150, x, 600-x ) 



3- 

x - 4-x 

4- (140625 -x 

#67: 

cl-e 

+ e -(c2-C0S(y34-x) + c3 

■ SIN 6/34 -x)) 

#68: 

SGC(5 , 

3 

-26, -150, x, 600-x ) 


#69: 

SGCCS, 

-x 

-26, -150, x, 20- e ) 


#70: 

Memory 

Exhausted 


#72: 

SGcfo, 

4, 0, x, ] 


#73: 

(. TAN (2-x) J 

LN(SIN(2-x)) C0S(2-x) • LN(TAN (x)) 

i 

- + c2-C0Sf2-x') + c3- SINl'2- 


140625 


LN(SIN(2-x)) CQ5(2-x) ■ LN(TAN(x)) 1 

#74: y (x) := + c2-C0S(2-x) + c3-SIN(2-x) + cl 


y rrr (x) + 4-y r (x) = C05(2-x) ■ 

#75 : SIN (2-x) 


5IN(2-x) J SIN(x) -COS(x) SIN(x) -COS(x) 








D-N-L#24 


J.L. Rodriguez & M.J. Fernandez: 3rd order ODEs 


p39 


#76: Simplify the difference of both right sides: 


#77: 


C0S(2-x) ■ 


SIN(2-x) 3 

SIN(2-x) J 


8-C0S(2-x) 


SIN(x) -COSCx) SIN(x) -COS(x) 
1 


3 TAN(2-x) 


#78 : 3 3 33 

SIN(2-x) SIN(x) -CGS(x) SIN(x) -CGS(x) 

#79: Trigonometry Collect 

#80: Tri'gpower := Si nes 

#81: £ 

#82: Result fram 1596: 

#83: C0S(2-x) ■ LN(TAN(x)) LN(SIN(2-x)) C0S(4-x) 


C0S(2-x) ■ 2- c2 


C0S(4-x) 
„ SIN(2-x) 


dx 


4 2 

#84: simplified: 

LN(SIN(2-x)) C0S(2-x)-LN(TAN(x)) 
#83: 


+ c2-C0S(2-x) + c3-SIN(2-x) + cl 


+ c3-SIN(2-x) + cl 

4 


3x 


9) y'" - 3 y" + 2 y' = 


1 + e 


3x 


a) Find a particular solution 

b) Find the solution with y( 1) = 2, y ’ (1) = 3, y" (1) = 4 


#87: SP 


3-x 


2-x 

#88: e 


-3, 2, 0, x, 

> 

1 + s j 

X 

LN(e +1) 3 


#89 : SPC 


-3, 2, 0, x, 


4 J 

3-x 

e 

x 

1 + e 


x f x 1 "] LN(e + 1) 

+ e -| LN(e ■+■ 1) - + 

2 J 2 


1. 2, 3, 4 


LN 

#90: 

X 

e +1 

2-x 

LN 

X 

e +1 

-2 

e 3 


v e + 1 j 

v, e + 1 j 

X 

+ e ■ 

LN 

K 

S +1 

2 

t 2 2 A j 

s. e + 1 J 


-1 1 
e -+- 2-e - — 

2 J 
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This file shows once more the power of a CAS in the hands of an experienced user and it 
raises once more the question about the necessity of so many drill examples in the text 
books for the future. And the future may have still begun. Josef 


p40 


Thomas Weth: 


A Lexicon of Curves 


(9) 


D-N-L#24 


Ebene Algebraische und 
Transzendente Kurven (9) 

Thomas Weth, Wurzburg, Germany 


Pascalsche Schnecken - Snails of Pascal 

„Wenn alle die Krafte, die bei der Entwickelung einer Pflanze mitwirken, mathematisch erkannt waren 
und ebenso der innere Mechanismus ihrer Organe, so wiirde man im stande sein die ganze Lebensentwi- 
ckelung durch Formeln darzustellen, insbesondere wiirde man die Gleichungen derjenigen Kurven erhal- 
ten konnen, welche den UmriG ihrer Blatter darstellen. Aber umgekehrt, wenn man auch diese Gleichun- 
gen kennte, wiirde man dennoch nicht das Leben jener Pflanze durch Formeln darstellen konnen; doch 
auch von diesem Ziele ist man noch weit entfemt indem man sich begniigen muG, die Blattumrisse 
durch Gleichungen darzustellen, die nicht exakt sondem nur in einfacher Weise angenahert diese wieder- 
geben“ (Foria, 1902, S.307). Bodo Habenicht versuchte ausgangs des 19. Jhdts. dem Geist der Zeit wis- 
senschaftlicher und technischer Hochstleistungen (Rontgenstrahlen, Eiffelturm) entsprechend auch die 
Botanik durch die Sprache der Mathematik zu erfassen (vgl. Die analytische Form der Blatter , Quedlin- 
burg, 1895). Eine der einfachsten Kurven, die sich in modifizierter Form unter den von Habenicht ange- 
gebenen „Blattkurven“ findet, ist die Pascalsche Schnecke, die auf Etienne Pascal, den Vater des bekann- 
ten Blaise Pascal, zuriickgeht. Wie sich herausstellt, konnen spezielle Pascalsche Schnecken zur Dreitei- 
lung des Winkels verwendet werden und spielen bei mechanischen Problemen ein Rolle, weswegen sie 
urspriinglich auch unter dem Namen "Sauveur's und de rHospital’s Zugbriicke" bekannt waren. Bereits 
vorher hatten sich auch Johann und Jacob Bernoulli mit dieser Kurve beschaftigt. 

"If all the forces which are contributing to the evolution of a plant would be recognized mathe- 
matically and also their internal mechanism of their organs then we were enabled to represent 
the whole process of its life by formulae. We specially would obtain the equations for the curves 
which form the contour of their leaves. But reversely, if we even knew the equations, we would 
yet be unable to express the plant’s life by formulae. This goal is very far as we have to be sat- 
isfied with equations which are reproducing the leaves’ contours only in a very approximative 
way. (Loria, 1902, p.307). At the end of the 19 th century Bodo v. Habenicht inspired by scientific 
and technical supreme achievements (X-rays, Eiffel tower) tried to describe botany by the 
mathematical language. ( The Analytical Form of Leaves, Quedlinburg, 1895). One of the sim- 
plest curves which can be found among Habenicht's “Leave Curves" is the Snail of Pascal, 
which goes back to Etienne Pascal, father of the well known Blaise Pascal. It appears that spe- 
cial Snails of Pascal can be used for the trisection of an angle and they play a role in several 
mechanic problems. For that reason they were called "Sauveur’s and de I'Hospital’s Draw- 
bridge". Before that time Jacob and Johann Bernoulli had dealt with this curve, too. 
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Konstruktion - Construction 


Genauso wie die Konchoiden des Nikomedes (vgl. 
Folge 5) lassen sich Pascalsche Schnecken kon- 
struieren - nur verwendet man nicht wie dort eine 
Gerade, sondem einen Kreis als Leitlinie. 
Kreiskonchoiden ergeben sich nach folgender 
Konstruktionsvorschrift: 

Gegeben sind ein Kreis mit Mittelpunkt M und 
Durchmesser b und ein Punkt O auf der Kreis- 
linie. Von einem Kreispunkt Q aus tragt man auf 
der Geraden OQ in die beiden moglichen Rich- 
tungen jeweils eine Strecke konstanter Lange a ab; 
die Endpunkte Pj und P 2 dieser Strecken sind dann 
Konchoidenpunkte zum gegebenen Kreis. 



Pascalsche Schnecken als Kreiskonchoiden far drei 
verschiedene Abstande a, erstellt mit dem TI 92. 


The snails can be constructed similar to the Conchoids of Nikomedes (Lexicon #5) - with a circle as di- 
rectrix instead of a line. 

Given is a circle with centre M and diameter b and a point O on the circle line. Take any other point Q on 
the circle line, draw the line OQ and find the two points P 1 and P 2 with QP 1 = QP 2 - a - const. Then P 1 
and P 2 are two points of the curve. 


Herleitung der Kurvengleichung - Dervivation of the equation of the curve 


Aus der Konstruktion ergibt sich fur Pi und P 2 : 

OPi = OQ + QPi bzw. r = b cos q> + a und 
OP 2 = OQ + QP 2 bzw. r = b cos <p— a. 

Lasst man fur r negative Werte zu, vereinfacht sich 
die Darstellung zur allgemeinen Polardarstellung 
Pascalscher Schnecken: 

r = b cos (p+a. 

This is the polar form. See the construction. 



Mit r = ^x 2 + y 2 berechnet man (mit 
DERIVE) sofort: 

a 2 (x 2 +y 2 ) - (x 2 +y 2 -b-x) 2 = 0 
Damit sind Pascalsche Schnecken symmetri- 
sche algebraische Kurven vierter Ordnung. 


Fur a = b erhalt man aus der Kurvenschar wie 


in der nebenstehenden Abbildung die Kardioi- 
de (vgl. Folge 7). 



We substitute for r =fx 2 + y 2 then we can see that these curves are symmetric algebraic 
curves of order 4. For a = b we obtain the Cardioid. Compare with Lexicon #7. 
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Weitere Erzeugungsweisen 

Es sei noch erwahnt, dass Pascalsche Schnecken 
sich auch durch andere Konstruktionen erhalten 
lassen: 

• als Rollkurve 

Rollt ein Kreis auBen auf einem festen Kreis ab, 
so beschreibt ein markierter Punkt auf dem rol- 
lenden Kreis eine Pascalsche Schnecke. 

• als Inversionskurve 

Bildet man einen Kegelschnitt durch eine Inver- 
sion an einem Kreis ab, dessen Mittelpunkt mit 
dem Brennpunkt des Kegelschnitts zusammen- 
fallt, so erhalt man Pascalsche Schnecken. 

• als Ortslinie merkwiirdiger Dreieckspunkte 
Betrachtet man zu einem gegebenen Kreis Seh- 
nendreiecke, deren eine Ecke A festliegt und bei 
dem der zugehorige Winkel a konstant ist, so ist 
der Ort der In- und Ankreismittelpunkte aller 
derartigen Dreiecke eine Pascalsche Schnecke. 

• als Ortslinie eines Winkelscheitels 

Bewegt sich ein konstanter Winkel so, dass 
seine Schenkel zwei feste Kreise standig beriih- 
ren, so beschreibt sein Scheitelpunkt eine Pas- 
calsche Schnecke. 


Other ways to obtain a Cardioid 

It is worth to be mentioned that we can ob- 
tain a Snail of Pascal in some other ways: 

• as Trochoid 

A point on a circle rolling outside on a fixed 
other circle describes a Snail of Pascal 

• as a Curve of Inversion 

If we map a conic by an inversion at a circle 
with its centre lying in a focal point of the 
conic, then we again receive a Snail of 
Pascal. 

• as locus of remarkable points of a tri- 
angle 

Observing triangles inscribed in a given 
circle with a fixed vertex A and the accom- 
panying constant angle a, then the locus of 
the centres of all the incircles and excircles 
results in a Snail of Pascal. 

• as locus of a vertex of an angle 

A constant angle moving with its sides 
touching two fixed circles gives a Snail of 
Pascal as the locus of its vertex. 


Dreiteilung des Winkels mit Pascalschen Schnecken - Angle trisection using S. of P. 

Gegeben sei eine spezielle Pascalsche Schnecke mit a = bl 2. 

Den zu drittelnden W in kel 1 co= ZSMP tragt man wie in nebenstehender Zeichnung an. P sei der Schnitt- 
punkt des einen Schenkels von <ymit der Pascalschen Schnecke. Dann gilt: 


Da das Dreieck QOM gleichschenklig ist: ZMOQ = ZMQO =: cp. 


Der Winkel ZQMT ist AuBenwinkel zum Dreieck QOM, hat also die GroBe 2<p. Im gleichschenkligen 
Dreieck PMQ gilt fur die Basiswinkel: ZQPM = ZQMP — — = 90° - y- • 


Fur den Winkel ZTMP gilt also: ZTMP = co + 90° und weiters 


ZTMP = ZQMT +ZQMP =2 (p + 


r o (P ^ 

90° -T- 

2y 


= 90° + -r"- 


Daraus folgt: a> = — bzw : 

2 3 2 


1 


Im folgenden wird nicht unterschieden zwischen Winkel und WinkelmaB. Die jeweilige Bedeutung ergibt sich aus dem 
Kontext. In the following we will not differ between angle and its measure. 
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Halbiert man also den Winkel ZMOQ = tp. so erhalt man ein Drittel des gegebenen Winkels co. 

Let's take a special Snail of Pascal 
with a = b/2. 

co = ZSMP is the angle to be trisected. 
P is the intersection point of one of its 
leg with the curve. The following can 
be deduced: 

As AQOM is isosceles we find ZMOQ 
= ZMQO = <p. ZQMT is an exterior 
angle of the triangle QOM, hence 
equals 2cp. APMQ is another isosceles 
triangle with ZQMP = ZQPM = = 90° 

_ £ 

2 ' 



ZTMP = co + 90° and ZTMP = ZQMT + ZQMP =2 cp + 


90° - 


P 

2 , 


= 90° + 


3 cp 


This leads obviously to co =— ,co =^~, and to — = — . 

2 2 3 2 

So if we bisect the angle ZMOQ = cpwe obtain a third of the given angle co. 


Anmerkung zum WorldWideWeb - Interesting WWW - pages 

Highly recommended: 

The MacTutor History of Mathematics archive 

http://turnbull.mcs.st-and.ac.uk/~historv/ 

Index of Biographies: http://turnbull.mcs.st-and.ac.uk/historv/Bioqlndex.html 

History Topic Index: http://turnbull.mcs.st-and.ac.uk/historv/lndexes/HistorvTopics.html 

Famous Curves Index: http://turnbull.mcs.st-and.ac.uk/historv/Curves/Curves.html 

Mathmeticians of the Day: http://turnbull.mcs.st-and.ac.uk/historv/Dav files/Now.html 

Additional Matrials: http://turnbull.mcs.st-and.ac.uk/historv/lndexes/Extras index.html 

Two screen shots of PLOMAT(house) from the TI-92 to plot matrices of points. 

(left side: Points discrete and large; right side: Points Connected) 


fFlTHOY Y FS Y FH Y FE^ Y FfiT p \ 

Zoom Trace ReGraph Math DrawK y \ 

fFiT^OY Y F3 

t f— Zoom Trace 

FH Y FEt Y FfiT Yf? *Y 1 

ReGraph Math DrawK y \ 

y 

□ 

□ 

y 

/ 

\ 


□ 





X 





TIFF: G RAD EXACT FUNC 

TIFF: G RAD EXACT FUNC 
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Two other approximations to square the circle 

Alfonso J. P oblation Saez, Valladolid, Spain 


3. A method from Poland 

This third method is due to the Polish Jesuit Reverend Adam Kochanski, in 1685. He was the first 
to use a steel spring in the suspension of a clock’s pendulum. 


Kochanski 's con- 
struction with a circle 
of radius r - 1. 


This reminded me 
on my school time - 
on the screen of a 
TI-92. Josef 


Start with a circle (centre = O and radius r). Proceed with a circle of radius OA and centre A (on the 
diameter BOA in order to find C. With centre in C, draw one more circle having the same radius and 
obtain D. Consider the segment OD and its intersection with the tangent to the initial circle passing 
through A. This gives you point E. Then F is on this tangent and verifies that EF = 3 * OA. The length 
BF is approximately OA*7i. What is here the value of 7i? 

4. Ramanujan's contribution 

Finally, I chose the construction of a great Hindu mathematician. Srinivasa Ramanujan. He gave 
us a lot of interesting formulae, one of which is implemented in DERIVE to approximate n. (see [5]). 


From a given circle, consider M, the midpoint 
of OA and T such that OT = j- OB. Then P is 

on the circumference such that TP is perpen- 
dicular to AB and Q is such that BQ = TP. Take 
S as the midpoint of AQ and D satisfies 
AD = AS. The segments TR, BQ and OS must 
be parallel. Now draw AC being tangent to the 
circumference and length equal to RS. Finally, 
BE = BM and X such that EX is parallel to CD. 
What is the relation between BX and re? 
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Probability Distributions 
Proof and Computations (1) 

Peter Mitic, Medstead, UK 


ABSTRACT 

The use of DERIVE to perform relevant algebraic computations in the context of 
probability distributions is discussed, with particular reference to some common 
distributions. Problems with using DERIVE in this way are noted, and some partial 
solutions are suggested. 

INTRODUCTION 

Certain components of courses in probability and elementary statistics can benefit from the algebraic 
manipulation facilities offered by DERIVE. The way in which the DERIVE interface is used (i.e. select 
an expression and operate on it) makes it particularly suitable for computations which involve sequen- 
tial operations. Consequently, the user can concentrate on wider aspects of problem solving and is 
free to appreciate the overall strategy used. In this paper we use DERIVE to prove certain results in 
probability theory and show how some shortcomings of DERIVE may be overcome. 

We consider how DERIVE may be used to obtain some standard results in probability theory, and 
assess its efficacy in doing so. The initial discussion centres on the ability of DERIVE to calculate 
means and variances, given some standard probability distributions. This involves summing series 
and evaluating integrals. It is an advantage, from a didactic point of view, to be able to do such com- 
putations directly. Looking up standard results is not meaningful unless it is accompanied by a good 
conceptual understanding. Routine computations can help to provide this. 

MEANS AND VARIANCES OF DISCRETE RANDOM VARIABLES 

We consider the cases of the binomial and geometric distributions, because the principal pedagogic 
and technical points are covered by these distributions. 

For a discrete random variable, X, defined in a domain, S, the definitions below may be used to calcu- 
late the mean and variance, //and a respectively, of X. 
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Comment of the Editor: I leave this contribution in its original form from 1996. Many func- 
tions which Johann developed have been implemented since then. You can find the 
DERIVE 6 - suitable file at the end of Titbits (9). titbits9_new.mth is the respective file. 

Titbits from Algebra and Number Theory (9) 

by Johann Wiesenbauer, Vienna 

Time is flying, isn’t it? It doesn’t seem so long ago that I wrote the first column in this series using 
DERIVE 2.56 then. In the meantime numerous new versions of DERIVE have come out taking ac- 
count of the various users’ wishes. The present culmination of this development is DERIVE for Win- 
dows 4.0 (called DfW in the following for short) which has become available these days. 

What I think about DfW? Well, despite a lot of apparent improvements such as the lightning-fast 
graphics, the new file management system or the enhanced scrolling capabilities, it isn’t exactly true 
that it made me feel all warm inside from the very beginning - you can’t teach an old dog new tricks, 
as the saying goes. What I found most irritating was the input window that pops up every few seconds 
- why on earth didn’t they use the main window for this purpose like every other CAS I know of ? -, 
but after minimizing it, i.e. increasing the screen resolution and maximizing the main window at the 
same time, I finally got used to it. Currently I am discovering new nice features of DfW every day, 
and it is no longer inconceivable that I end up in saying one day: “I like it”. 

One of those nice features I was talking about is the option to include the main variable as a third pa- 
rameter in the functions REMAINDER and QUOTIENT. (Up to now, the main variable could only be 
adjusted by using “Manage Order” from the main menu.) To the best of my knowledge, this was sug- 
gested by Eugenio Roanes-Lozano, who had to apply in [3] a lot of trickery to circumvent this restric- 
tion when defining his function PREM(p,q,v): 

POLY_COEFF(ax # n): = l/n!*LIM(DIF(ax # n) f x # 0) 

POLY_DEGREE_AUX(u / x,n): = IF(u=0 / n / POLY_DEGREE_AUX(DIF(u / x) / x / n+ 1 ), 
POLY_DEGREE_AUX(DIF(u,x) / x / n+ 1 )) 

POLY_DEGREE(u / x): = POLY_DEGREE_AUX(u / x / -l ) 

LCOEFF(p / v): = POLY_COEFF(p / v / POLY_DEGREE(p / v)) 

arb_var:=(RHS(SOLVE(0=0 / x))) SUB 1 

SUBST(f,x,a): = LIM(f,x,a) 

SWAP_AUX(f / x / t / y):=SUBST(SUBST(SUBST(f / x / t) / y / x)4 / y) 

SWAP(f,x,y) : = SWAP_AUX(f / x / arb_var / y) 

REMAINDER_V_AUX(p / q / v / vO):=SWAP(REMAINDER(SWAP(p / v / vO) / 

SWAP(q,v,vO)),v,vO) 

REMAINDER_V(p / q / v): = REMAINDER_V_AUX(p / q / v / (VARIABLES([p / q])) SUB 1 ) 

MUL_FAC(p / q / v): = LCOEFF(q,v) ~ (1 +POLY_DEGREE(p / v)-POLY_DEGREE(q / v)) 

PREM(p / q / v): = REMAINDER_V(MUL_FAC(p / q / v)*p / q / v) 
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In DfW the same function can be defined without any auxiliary functions (!) as follows: 

PREM(p,q,v):=REMAINDER(LIM((TERMS(EXPAND(q,v))) SUB l,v,l)~ 

(1 + LIM(v*DIF(f_,v)/f_,f_,(TERMS(EXPAND(p,v))) SUB 1)- 
LIM(v*DIF(f_,v)/f_,f_,(TERMS(EXPAND(q,v))) SUB l))*p,q,v) 

On close inspection of this code you will find out that I used the following more efficient implementa- 
tions 

POLY_DEGREE(u < x): = UM(x*DIF(f_,x)/f_,f_,(TERMS(EXPAND(u,x))) SUB 1 
LCOEFF(u,x):= LIM((TERMS(EXPAND(u,x))) SUB l,x,l) 

Though I didn’t need POLY_COEFF(u,x,n) I would like to point out the following interesting alter- 
native to the definition given above (cp. DNL #21, p37): 

POLY_COEFF(u,x,n):=QUOTIENT(REMAINDER(u,x^ (n+1 ),x),x^ n,x) 

Furthermore, for those with older versions of DERIVE the following implementation of PREM might 
still be of interest (if @ is not available in your DERIVE version then use any other unused variable 
instead of it!): 

SWAP(f,x,y): = ITERATE(LIM(f_,@,y),f_,LIM(f,[x,y],[@,x]),l) 
REM(p,q,v):MTERATE(SWAP(REMAINDER(SWAP(p,v,vJ,SWAP(q,v,vJ),v,v_),v_, 

(VARIABLES([p,q])) SUB 1,1) 

PREM(p,q,v):=REM(LIM((TERMS(EXPAND(q,v))) SUB l,v,ip(l +LIM(v*DIF(f_,v)/f_,f_, 
(TERMS(EXPAND(p,v))) SUB 1 )-UM(v*DIF(f_,v)/f_,f_, 

(TERMS(EXPAND(q,v))) SUB 1 ))*p,q,v) 

You might wonder by now what this mysterious PREM-routine is all about. Sorry, I should have told 
you that earlier! Above all it is a very useful tool to solve a system of polynomial equations by trans- 
forming it into a solvable system of equations without “losing” any solutions of the original system. 
To be more precise, if / and g are polynomials in x,,...,x n , where we may assume that 
deg(/,Xj) > deg(g,x,) w.l.o.g., then either deg(g,Xj)) = 0 or we can shift to polynomials q, r in 
x i ,...,x n , which are defined by the following equation (reflecting the polynomial division of m -/by 
&): 

m-f = q-g + r with deg(r,Xj )<deg(/,Xj). 

Here m is a polynomial inx 1? ...,x n , that is chosen in such a way that fractions are avoided (see the 
formula MUL_FAC(p,q,v) above for its exact definition) and the pseudo-remainder r is the output of 
PREM(f,g,xl). Again either deg(r,x 1 ) = 0 or we can repeat this procedure with g and r until we finally 
arrive at polynomials /', g' , s.t. g f doesn’t contain x\ any longer, and it is clear from the construction 
that every solution of “/= g = 0 is also a solution of f r = g r = 0. 

Given a system of polynomial equations this algorithm can be used to achieve a “triangulated system” 
in a way very similar to GauB’ elimination. I will try to give you an idea how this works by means of a 
simple example, which was also used by A.Perotti in [ 2 ]. We consider the system of polynomial 
equations 
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x-y- 1-z 
x-z — l — y 
y-z- l-x 

By applying PREM three times 


( 1 ) 


PREM(x*y-l +z,x*z-l +y,z)=x' v 2*y-x-y+l (2a) 

PREM(x*y-l +z,y*z-l +x,z)=x*(y^ 2-1 )-y+l (2b) 

PREM(x*(y' v 2-l)-y+ 1 ,x^ 2*y-x-y+ 1 ,y)=-x~ 5+x^ 4+2*x~ 3-3*x~ 2+x (3) 

we get the triangulated systems made up of (l),(2a),(3) or (l),(2b),(3), respectively. Solving (3) w.r.t. 
x yields 

SOLVE(-x ^ 5 +x ^ 4 + 2*x ^ 3-3*x ^ 2 +x,x) = 

[x=0,x= 1 ,x=SQRT(5)/2-l/2,x=-SQRT(5)/2-l/2] 

Using (2a) and (2b) we get the corresponding values of y: 

SOLVE(LIM(x ~ 2*y-x-y+ 1 ,x,0),y)=[y= 1 ] 

SOLVE(LIM(x^ 2 y - x - y + 1 , x, 1 ), y) = [y = @2] 

SOLVE(LIM(x*(y ~ 2-1 )-y+ 1 ,x,l ),y)=[y=0,y= 1 ] 

SOLVE(LIM(x^ 2*y-x-y+l ,x,SQRT(5)/2-l/2),y)=[y=SQRT(5)/2-l/2] 

SOLVE(LIM(x^ 2*y-x-y+l ,x,-SQRT(5)/2-l/2),y)=[y=-SQRT(5)/2-l/2] 

Now for reasons of symmetry y = z must hold except for x = 1 , where we only know that y,z s{0,l}. 
By using (1) to decide this ambiguous case we finally arrive at the five triples 




All these triples are really solutions of the original system. Note that we have actually used 
(1), (2a), (2b), (3) to compute them. It is easy to see that both triangulated systems (l),(2a),(3) and 
(l),(2b),(3) on their own have more solutions than the original system. (Check it!). 


One could solve this system of polynomial equations also by means of my SOLVE2-function from my 
last ‘Titbits’ but it is very tricky and cannot be generalized. I hate to admit it, but Eugenio’s PREM- 
function is streets ahead of my RED-function (the kernel of SOLVE2), when it comes to solving sys- 
tems of polynomial equations. (It must be said though that it was invented for quite a different pur- 
pose, cf. DNL #20, p38.) If you are interested in learning more about applications of the PREM- 
function I can recommend to you reading Eugenio’s wonderful paper [3] where he used it extensively 
to prove theorems of elementary geometry by the so-called Wu’s method. At any rate I don’t want to 
leave this topic without thanking him for his cooperation and many fruitful discussions via email. 
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Back to DfW. What else could be done to enhance it? Well, a lot of little things come to my mind, e.g. 
that there is no calculation time shown after the instant simplification of expressions by “=” or the fact 
that there is no Cartesian product for lists and sets, and what not. But from a number theoretic point of 
view there is one thing that is on the top of my wish list and once more it concerns the notorious pow- 
ermod- function: It should also work for negative exponents! 

At present if you ask e.g. for the inverse of 3 mod 5, i.e. mod(3 A (-l),5), you will get the answer 
mod(3 A (-l),5)=l/3, 

which is not exactly what you want. It could be programmed along the following lines (cp. DNL #17, 
p35): 

POWERMOD(an / m): = IF(n>=0 / MOD(a /v n / m) / IF(GCD(am)=l / MOD(MOD((ITERATE( 
IF(MOD(a_ / b_)=0 / [a_ / b_ / c_ / dJ / [b_ / MOD(a_ / b_) / d_ / c_-FLOOR(a_ / b_)^d_]) / 
[a_ / b_ / c_ / dJ / [a / m / l / 0])) SUB 4,m) A (-n),m))) 

If the exponent n is negative and gcd(a,m)>l, the output will be a question mark as, of course, power- 
mod^, n,m) is undefined in this case. 

I will use this function to implement another very important formula which is based on the so-called 
Chinese Remainder Theorem (CRT): 

If m 1? m 2 ,...,m r is a sequence of pairwise coprime integers, then for arbitrary integers a x ,a 2 ...,a r the 
system of congruencies 

x = a x mod m x , 
x = a 2 mod m 2 , 


x = a mod m , 

has a solution, which is unique mod .m r . 

Fortunately, there exists also an explicit formula for the solution which can be written in DERIVE - 
notation as follows (cf. [1], p37) 

CRT(a,m): = ITERATE(MOD(a*VECTOR(p_/m_*POWERMOD(p_/nn_ # -l ,m_) / nn_ / nn) / p_) / 
p_,PRODUCT(m_,m_,m),l) 

Here a denotes the vector [a x ,a 2 ...,a r \ and m the vector [m 1 ,m 2 ,...,m r ] . Again a question mark will be 
output if the numbers m 1 ,m 2 ,...,m r are not pairwise coprime. Number theory abounds with applica- 
tions of this important routine, but space is mnning out. Therefore, bye for now! (email: 
i . wiesenbauer@tuwien.ac.at) 
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As I don’t want to print all the Tl-programs and -functions which I have received for the DNL's Tl - 
Corner - many thanks to you all -I include the files in a special Tl-directory on the Diskette of the Year. 
I am sure that not all Tl users have Tl-Graph LINK at their disposal so they will not be able to use the 
programs. I tried to find another way: I include a document (in Word for WINDOWS 6 format) contain- 
ing many programs and functions in readable form. (User Forum, Terence Etchells) You will find W. 
Propper’s wonderful DIRA - for discussing curves - in a German and in an English version (DNL#23). I 
also have the pleasure to enclose D. Stoutemyer’s Program Library together with a printable docu- 
mentation. Many thanks David. DIRA and the library are not on a Word document because they are 
too big. At last I want to point out that there is lot of useful Tl-stuffto be downloaded from the Internet. 
(Web pages of Tl, SWH, ....) 


At last I'm glad to announce a lot of Tl-contributions for the next year. Most of them will be interesting 
for DERIVE Users, too. (So I think that E.Laughbaum's contribution for the Tl on the next pages can 
easily be transferred on PC-DERIVE, see page 47.) There will be among others: 

Dynamic System, Permutations, Numbers of Bernoulli, Continued Fractions, R.Schorn; Mortgage Ta- 
bles, Tania Roller; The Simplex Method on the Tl, Bruce Chaffee; Diophantine Equations & The 
Gauss-Seidel Method, L.Tortosa & J.Santacruz; Endpoint vs. Interior Extrema, White & Leinbach; 
Investigations on GCD and LCM with the Tl, Griebel; First experiences with a Tl in the class room, 
Bohm.... 


SOLSYST - A Function to solve Simultaneous Equations 

BY W. PROPPER, NURNBERG, GERMANY 


Sol syst (sys , var) 

Func 

Local i , j , n , mtx , dtm 
rowDim(sys)-»-n 
newMat ( n , n+1 ) -*mtx 
For i , 1 , n 
For j , 1 , n 

^(1 eft (sys [i ,1 ] ) ,var[j ,1 ] )-nntx[i , j ] 
ri ght (sys [i , 1 ] )->-mtx [i , n+1 ] 

EndFor 

EndFor 

det (subMat (mtx , 1 ,1 ,n,n))+dtm 

If getType(dtm)="NUM" and dtm*0 or getType(dtm)="EXPR" Then 
rref (mtx)-»-mtx 

Return subMat (mtx , 1 , n+1 , n , n+1 ) 

Else 

Return "Nicht eindeutig losbar" 

Endlf 

EndFunc 
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With this nice function you can solve systems of linear equations in a DERIVE like way. Non 
German speaking users have to change one string in the third line from the bottom: "Nicht 
eindeutig losbar" to the equivalent of "No unique solution". 
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Sums of Absolute Value Functions: 
An Application 

Edward D. Laughbaum, Columbus, Ohio 
elaughba@math.ohio-state.edu 


Consider the following: 

a) The level of the drug Imipramine in the blood of a patient rises at a constant rate (for example 60 
nanograms per week) until the patient is at the prescribed level and then the rate of change re- 
mains at 0% until the patient is taken off the drug at a constant rate (for example 90 nanograms 
per week). 

b) For users of electricity who also use a heat pump, many power companies charge, for example, 
$0.08 per kWh used for the first 1000 kWh and $0.05 for the next 1000 and finally, $0.03 for any 
consumption over 2000 kWh. 

c) Sliding scale commissions are business examples of multi-constant rates of change and are used 
to provide incentive for higher employee productivity. For example, a business may pay the sales 
staff 3% commission on sales from $0 to $10,000, 5% on sales from $10,001 to $15,000, and 8% 
commission on sales $15,001 and over. 

d) The pay scale for piece-work by a telemarketing company is $0.35 for the first 300 calls in the 
week, $0.42 for the next 200 calls, and $0.65 for any call over 500 calls. 

e) A new long-distance phone company has the following rate schedule: 

120 per minute for the first 15 minutes 
90 per minute for the next 10 minutes, and 
60 per minute for any time over 25 minutes. 

f) The 1994 (or any year) United States federal income tax form 1040 schedule (for single filers) 
had rates of taxation of 15% on the first $22,750 of taxable income, 28% on the next $32,350, 
31% on the next $59,900 of taxable income, etc. That is, the taxable income brackets are at 
$22,750, $55,100, and $1 15,000. (There were two more brackets that will be ignored for the sake 
of brevity.) 

g) A commercial airline flight from Reno, Nevada to St. Louis, Missouri ascends at a constant rate 
(after initial take-off) of 1 100 feet per minute for 30 minutes (until it reaches a cruising altitude of 
33,000 feet). It then levels off until it is 60 minutes into the flight. At 60 minutes into the flight, it 
has burned off enough fuel to ascend to 37,000 feet at a rate of 400 feet per minute. This takes 10 
minutes. The plane remains at 37,000 feet until 175 minutes of flight time when it descends at a 
rate of 1233 feet per minute and then lands in St. Louis. 
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The Algorithm for Finding the Model. 

The algorithm assumes a knowledge of behavior of sums of absolute value functions (Laughbaum, 
1996). 

Recognize that the structure of the model for all of the above situations is 
M = a\x + e x | + b\x + e 2 j + ... d\x + e n | + / , where there are n corners. 

1 . Find the comers (e\, e 2 , . . .) [corners are normally known] 

2. Simplify the model for each rate interval using the TI-92 

3. Set the coefficients of x equal to each rate of change. [rates of change are normally known] 

4. Solve the system. 

5. Find / using a geometric transformation of a vertical shift. 

An Example: 

The on-board computer on a commercial airline flight from Reno, Nevada to St. 
Louis, Missouri commands the plane to ascend at a constant rate (after initial 
take-off) of 1100 feet per minute for 30 minutes (until it reaches a cruising alti- 
tude of 33,000 feet). It then levels off until it is 60 minutes into the flight. At 60 
minutes into the flight, it has burned off enough fuel to ascend to 37,000 feet at 
a rate of 400 feet per minute. This takes 10 minutes. The plane remains at 
37,000 feet until 175 minutes of flight time when it descends at a rate of 1233 
feet per minute and then lands in St. Louis. 

1 . Recognize that the stmcture of the model is 

h(x ) = a\x + e x | + b\x + e 2 \ + c\x + e 2 1 + d\x + e A | + e\x + e 5 1 + / , where there are 5 comers. 

2. Find the corner parameters (el, e2, . . .). They are given as 0, -30, -60. -70, and -175. Thus, the 
initial model is h(x) = a\x\ + b\x - 30| + c\x - 60| + d\x - 70| + e\x - 175| + / , where x is time in 
minutes. 

3. Simplify the model for each rate interval using the TI-92 

Conditions: 

x > 0 andx < 30 

x > 30 andx < 60 

x > 60 andx < 70 

x > 70 andx < 175 
x> 175 
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4. Set the coefficients of x equal to each rate of change; that is, each piece of the model simplifies to 
a linear function as shown on the TI-92 screens. Thus, the coefficient of x is the rate of change for 
each piece. 

a-b-c-d -e= 1100 
a +b - c- d - e - 0 
a+b + c-d — e — 400 
a+b+c+d -e = 0 
u-\-b-\-c-\-d-\-6 = —1233 

5. Solve the system. 




The model for the height control function is 


h(x) = -— |x| — 550|x — 30| + 200 |x - 60| - 200 |x - 70| — |x-17 

6. Find /using a geometric transformation of a vertical shift. 

When x is 0, for example, the function is -126387.5. Add 126387.5 as the parameter /. 



The final model needed by the computer is 


h(x) = x -550x-30| + 200|x - 60| - 200|x - 70| jc — 1751 + 1 26387.5 . 


References 

Laughbaum, E. D. (1996). Modeling data exhibiting multi-constant rates of change. 

The AMATYC Review, 17(2), 27-34. 

(You can find the equivalent DERIVE application as SUM_ABS.MTH. Josef) 










